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ABSTRACT 


i 


Problems related to an unmanned exploration of the planet Mars 
by means of an autonomous roving planetary vehicle have been investigated. 
These problems include: design, construction and evaluation of the vehicle 
itself and its control and operating systems. More specifically, vehicle 
configuration, dynamics, control, propulsion, hazard detection systems, 
terrain sensing and modeling, obstacle detection concepts, path selection, 
decision-making systems, and chemical a* alyses of samples have been studied. 

Although earlier work along these lines had been focused on a 
vehicle which would house, transport and service scientific equipment for 
on-site analysis, more recent work has been directed to a vehicle capable 
of gathering specimens and data for an Augmented Viking Mission or to pro- 
vide the basis for a Sample Return Mission. 

The following tasks have been under active study: design, con- 
struction and testing of a dynamic vehicle, (weight exclusive of scientif- 
ic payload of 150-200 pounds); mathematical modeling and experimental 
measurement of the vehicle dynamics; propulsion and deployment systems; 
electro-mechanical control systems; telemetry systems to permit autonomous 
off-board control; collapsibility concepts and hardware; design, construc- 
tion and evaluation of double-hoop toroidal wheels; design and evaluation 
of a triangulation-based hazard detection system; obstacle perception by 
edge enhancement techniques; terrain modeling based on gradients; obstacle 
classification; path selection systems simulations and evaluations; gas 
chromatograph systems design concepts; experimental gas chromatograph sepa- 
ration measurements; and gas chromatograph model developments. 

These tasks and other subtasks are defined in detail and the 
progress which has been achieved during the period July 1, 1974 to 
June 30, 1975 is summarized. 

The educational impact of this project for the year has been 
assessed. A detailed report on the education implications of this NASA 
project at Rensselaer over the past eight years is included as Appendix A. 
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. Analysis and Design of a Capsule landing System 
and Surface Vehicle Control System for Mars Exploration 


1. Introduction 


Current national goals in space exploration include a detailed study and 
examination of the planet Mars. The effectiveness of Mars exploration missions 
would be enhanced according to the extent to which the investigative devices which 
are landed are mobile, to the range of their mobility, and to the ability to control 
their motion. In order to achieve basic mission objectives, and beyond that, to 
maximize the return on the commitment of resources to the mission, formidable techni- 
cal problems must be resolved. The major factor contributing to these problems is 
the round trip communications time delay between martian and earth control stations 
which varies from a minimum of about 9 minutes to a maximum of 25 minutes. This time 
delay imposes stringent requirements on the vehicle, on its control and communication 
systems and on those systems included on board to make the scientific measurements, 
in terms of their ability to function autonomously. These systems must be able to 
operate with a high degree of reliability and must be capable of calling for earth 
control under appropriate circumstances. 

A number of imports’. t problems originating with these factors and relating 
directly the basic mission objectives of an unmanned exploration of Mars have been 
and are currently being investigated by a faculty-student project team at Rensselaer 
Polytechnic Institute with the support of NASA Grant NCL 33-018-091. 

In addition to the technical objectives, the project is intended to involve 
students in meaningful design and problem-solving activities and by so doing to en- 
hance their formal education in the area of design. 

This progress report describes the tasks which have been undertaken and 
documents the progress which has been achieved in the interval July 1, 1974 to 
June 30, 1975 in addition to appraising the educational implications of the project. 


II. Overview of the Project 


For a number of years, the thrust of the project was derived from the 
objective of autonomous roving laboratory capable of a range of at least 1000 kilo- 
meters and a minimum operational life of one year. More specif ically .Rensselaer 's 
goal was to conceive, design and evaluate a roving vehicle concept together with 
the necessary hazard detection and control systems which could house, protect, 
service and transport the scientific instruments required to gain the desired bio- 
logical, chemical and geological data. 

However, the Viking project experience has cast considerable doubt on the 
feasibility of a sophisticated roving laboratory for purposes of a thorough unmanned 
exploration of Mars. This is not to say that the Viking project will not meet the 
prime mission objectives but rather to suggest that the far higher level goals anti- 
cipated for subsequent rover missions may well be met far better by alternative 
strategies. 

The basic problem orginates with the extreme difficulty of miniaturizing 
scientific instruments and apparatus and of building in the minimal level of adapt- 
ability required to deal with unforeseen or unanticipated situations. 
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Two strategics which appear to have significant potential have been identi- 
fied: (a) a series of Viking-type landers each of which would contain specialized 

equipment for a particular set of scientific studies but none of which would attempt 
to house all of the equipment required for a thorough exploration, and (b) sample 
return missions which would bring back specimens to earth for detailed analysis in 
earth laboratories. At this time, the preference of the scientific community appears 
to be a sample return mission which has the advantage of using the very superior 
capabilities of earth-based laboratories for martian samples analysis. 

However, for either of these strategies to be viable, it has been concluded 
that some kind of a roving vehicle capable of gathering samples and returning them to 
the stationary lander will be required. Although the minimum capabilities of such a 
sample gathering rover have not been established an estimate of the range can be 
made. The size of the two sigma landing site footprint is of the order of 100 x 200 
kilometers. Furthermore the areas considered attractive by scientists after analysis 
of the Mariner Orbitor Missions are too dangerous for landing either of a Viking- type 
of lander or a sample return vehicle. Thus it would appear that a sample-gathering 
rover capable of a 300 kilometer radial displacement would be required. Depending on 
the nature of the terrain between the landing site and the sample procurement loca- 
tions, the displacement distance between these two, the capabilities of the rover 
and the effectiveness of hazard detection and path selection systems, the total dis- 
tance to be traversed could be as high as 1500 to 2500 kilometers. 

Given this overall situation, the thrust of the project was redirected in 
January of 1975 with the concurrence of the Office of Planetary Studies and the Jet 
Propulsion Lab towards the conception, design and evaluation of a mini-rover capable 
of augmenting Viking-type lander missions or enabling a sample return mission. The 
redirection of the project is more a matter of emphasis than of substance since the 
basic goals remain to develop a vehicle of exceptional mobility and reliability 
together with the hazard detection and path selection systems required for autonomous 
roving. A sample gathering mission imposes stringent requirements on a vehicle in 
terms of its ability to negotiate hazardous terrain, its maneuverability in avoiding 
truly hazardous situations, its hazard detection systems and its path selection 
decision making system. 

Since January 1975 the major emphasis has been directed to tasks relating 
directly to a vehicle designed for a sample and data gathering mission. Beyond 
this emphasis, work has continued in the area of gas chromatographic systems analysis 
and design. The following major tasks have been under study: 

A, Vehicle Configuration, Control, Dynamics, Systems and Propulsion 
The objective of this task is to conceive, design and evaluate 

a roving vehicle for data and sample gathering during an unmanned 
exploration of Mars. Included are the aspects of vehicle con- 
figuration, collapsibility for launch configuration, deployment, 
dynamics, and motion and attitude control. One goal is to 
develop a vehicle concept capable of operating within the con- 
straints of a Mars mission involving an augmented Viking or a 
sample return vehicle and of adapting to changes in the mission 
as a result of the information gathered during the exploration. 

The second goal is to construct a rover which can serve as a test 
bed for the evaluation of hazard detection and avoidance systems. 

B, Hazard Detection Systems , The objective of this task is to 
conceive, develop, construct and evaluate a short range (0,5 to 
2.5 meter) hazard detection system to provide data for an on- 
board path selection decision making system. 
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C. Path Selec tio n Systems . The objective of this task is the 
simulation and evaluation of path selection systems with parti- 
cular emphasis on the development of path selection algorithms. 

This task is also responsible for providing guidance for the 
design of hazard detection hardware, 

D. Obstacle Detection and Terrain Modeling , This task is concerned 
with the problem oi detecting obstacles and modeling terrain for 
the purpose of path selection. 

E. Chemical Analysis of Specimens . A major objective of martian surface 
exploration will be to obtain chemical, biochemical or biological 
information. Many experiments proposed for the mission require a 
general duty, gas chromatograph-mass spectrometer system for chemical 
analysis. The objective of this task is to generate fundamental 
data and concepts required to optimize this chemical analysis system. 

The project milestones for the past year and as planned for the period 
July 1, 1975 to June 30, 1976 are shown in Figure 1. 


III. Summary of Results 

Task A, Vehicle Conf iguration, Control, Dynamics, Systems and Propulsion 

This broad task has been subdivided into the following subtasks: experi 

mental and analytical studies of vehicle dynamics; detailed design of structural 
components ; design of deployment hardware; motor drive and mechanical transmission 
details; development of a deployment procedure; detailed wheel design and construe 
tion; on-board controls for steering, motor drive for levelizing, deployment and 
emergency maneuvers; radio link for remote operation; and remote control station 
and transmitter design. 

A. 1. Vehicle Dynamics 

On the basis of observed dynamical behavior of the vehicle during 
extensive field testing, a review of dampening alternatives was 
undertaken. A two-dimensional, two-degree of freedom dynamic 
model was used to select parameters for rotational dampening of 
the rear struts. The vehicle dynamical response has been much 
improved by the addition of rotational dampers. 

A. 2. Structural Aspects Including Propulsion 

Following extensive field testing during the summer and early 
fall of 1974, the vehicle was disassembled and the components 
were examined for wear and failure indications. It was determined 
that design changes in the front end and torsion bar assembly were 
desirable. These changes have been implemented. Sensors have also 
been mounted such that the orientation of all major structural 
members will be known relative to the payload box. The data provided 
by these sensors in concert with the attitude gyro data will provide 
unamblgious vehicle state parameters to the remote controller. A 
laboratory obstacle course has been constructed and can be used to 
evaluate the vehicle's ability to negotiate special terrain features 
and the validity of the structural component position sensors. 
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Early hopes that rear wheel drive would suffice were unsupported 
by field testing* Analysis of the vehicle mobility indicated 
that a substantial advantage could be gained with front wheel 
drive. Front wheel drive was added and its advantages in increas- 
ing the vehicle’s mobility in difficult terrains was demonstrated* 

A. 3. Wheel Design and Testing 

The toroidal wheel, which had been redesigned during the previous 
year to increase lateral stiffness was further modified to de- 
crease the footprint pressure and to increase its traction relative 
to irregular surfaces. A double-hoop toroidal wheel meeting de- 
sign objective?^ was constructed but has not yet been field tested 
on the vehicle. 

A . 4 . Telemetry Systems and Vehicle Electronics 

Three tasks were pursued during the past year. Vehicle electronic 
systems were updated and modified to achieve desired operational 
characteristics. . control system for four-wheel drive was de- 
signed and imp lemon’ d. The most significant effort was in the 
design and construction ot a data link to be interfaced ultimately 
to a computer for autonomous vehicle roving. 

A 50 Kbit/sec telemetry link from the vehicle to the controller 
receiver hes been constructed, installed on the vehicle and tested. 
At the present time, this link transmits the strut sensor, steering 
angle, motor speed, latch and attitude gyro data to the off-board 
receiver for display. This link will be able to handle the hazard 
detection system data required for hazard avoidance decisions by a 
computer. 

The original manual control box was modified to permit more effec- 
tive remote control by an observer. This system is required for 
the extensive field testing program planned prior to autonomous 
roving studies. 

Some of t he software required for autonomous roving has been 
developed, CRT display of the front, side and top view of the 
vehicle including its pitch and roll is one of the software objec- 
tives. 


Task B, Hazard Detection Systems 

With the new objective of achieving a minimal level of autonomous 
roving by April 1976, the effort was directed to a short range 
(0.5 - 3.0 meter) hazard detection system. A concept employing a 
solid state laser and a solid state detector and using triangula- 
tion was constructed and demonstrated to be feasible as a detector 
of hazards to be avoided. 


\ 
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Task C. Path Selection Systems Simulation, Development and Evaluation 

Two potential path selection systems received in-depth simulate 
and evaluation. The Mid-Range System which employs two or to. e 
body-fixed scanners has indicated real potential for the do vCt<on 
of positive obstacles. The Short Range System (1.0 to 2.0 meters 
range) has demonstrated an a>ility to detect craters providing, 
therefore means for avoiding negative obstacles. Although it is 
not now planned to implement either of these concepts directly, 
the knowledge gained will aid the implementation of the hazard 
detection concept developed under Tark B. Modeling of the latter 
concept has been initiated along wit'll more detailed terrain simula- 
tion capabilities so as to per t a full scale evaluation of this 
concept as a potential hazard detection system. 


Task D. Obstacle Detection and Terrain Modeling 

This task deals with the problems of obstacle detection, terrain 
modeling, and scanning schemes for data acquisition. 

D.l.a. Obstacle Detection and Picture Enhancing 

The Four Directional Ratio Algorithm developed previously was 
able to detect only the tops of boulders and the leading edges of 
craters. A new algorithm making use of a Kalman Filter has been 
developed which not only detects the tops of boulders and leading 
edges of craters but also the bottom edges of boulders and the 
far edges of craters. 

D.l.b. Minicomputer Simulation of Obstacle Edge Detection Systems 

A computer program has been developed for use with a graphics 
terminal minicomputer to simulate and display the output of obstacle 
edge detection systems. The program allows the investigator to 
perceive visually the effects of data density and accuracy, obstacle 
shape and location on the identification of obstacles by edge detec- 
tion and to adjust parameters readily. 


D.2. 


Obstacle and Terrain Feature Classification 


The problem of identifying the minimum number of terrain parameters 
required to deduce the passability or impassability of terrains has 
been examined. It is concluded that at least four parameters and 
three inequality conditions which must be met simultaneously are 
required. These results will be valuable input to the groups working 
on hazard detection hardware alternative and path selection algorithms. 


D.3. Terrain Modelin g 

A two-dimens ioni* » spline function procedure has been developed for 
the modeling on terrain using innge and pointing angle data such as 
may be provided by a laser scanning system. The modeling procedure 
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fits low order two dimensional polynomials to two adjacent terrain 
elements and provides a basis for minimizing discontinuities in 
height and slope at the common boundary. Once the first two terrain 
elements are modeled, a third adjacent element can be modeled with- 
out change in the parameters describing the first two elements. This 
process can now be extended to any number of terrain elements and 
ultimately yields gradient information to be used in a path selection 
algorithm. This terrain modeling procedure is intended to complement 
ths obstacle detection concepts, Ta3k C.l.a. 

! 

D.4.«a. Measurement Scanning Concepts 

Three alternative concepts have been considered: (a) Terrain and 

Zoom Scan, (b) Obstacle Scan and (c) Corridor Scan. The number of 
data points required for each method have been determined and the 
length of time required assuming onr millisecond per measurement has 
been estimated. The relationship between vohiclr velocity, maximum 
scanning range and time for data acquisition has been studied. It 
appears that the time for data acquisition can vary from about 6 to 
15 seconds assuming a one millisecond data acquisition interval. At 
high vehicle speeds of the order of 1 meter/second, a scan range of 
30 or so meters may be required for a timely detection of hazards. 

At low speeds of 0,1 me ter /second , the scan range and concomitant 
data requirements can be reduced significantly without exposing the 
vehicle to hazardous situations. 
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D.4.b. Laser Rangefinder Errors 

The implications of spread inp of the laser bean and of the resulting 
reflections being sensed by a detector as possessing different times 
of flight and therefore different ranges have been investigated. It 
is shown that for a laser mounted 3 meters above the ground that the 
optimum laser aperture at the 30 meter range is 0.78 centimeters. 

For this set of parameters, the range error duo to beam spreading on 
a horizontal terrain is 6.27 cm. Although this error is predictable 
for a level terrain and could conceivably be corrected, the fact that 
a rover would encounter terrains of all possible orientation excludes 
this simple correction. 


Task E. Chemical Analysis of Specimens 

E. 1. Chromatograph System Characteristics 

These studies are intended to provide chromatograph design techniques 
and guides using simulation models currently under development. 
Specific studies include: (1) development of a generalized simulation 
proce J urc; (2) determination of the effect of design parameters upon 
the degree of separation; and (3) characterization of realistic sample 
injection pulses. The simulation procedure treats a maximum of five 
components anu uses actual sample sizes, carrier gas flow rates, and 
real time so results can be related to real systems. A method for 
quickly predicting degree of separation was developed by usiug a 
Cuussian-typc approximation to the chromatograms. Input pulse data 
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were obtained on the experimental test facility. First and 
second time moments were correlated to provide a means of predict- 
ing gross characteristics of such pulses. 

E.2. Chromatograph Simulation Development 

Models developed previously to explain nonlinear composition- 
dependent effects appear too complicated to be solved efficiently 
by numerical methods. A simpler model involving the nonlinear 
Langmuir isotherm was derived and was solved numerically by the 
Crank-Nicolscn algorithm. This model represented well a chromatogram 
of n-hept on Cnrotnosorb C- \02 previously simulated poorly by the 
linear model because of much tailing in the data. The numerical pro- 
cedure was inefficient, requiring 25 to 60 minutes of IBM 360/67 time 
to fit a set of chromatographic data, so future work will involve 
development of a more practical numerical procedure. 


IV. Educational Implications 

As noted earlier, one of the objectives of the project is to involve 
engineering students in meaningful design and problem-solving activities enhancing, 
therefore, their overall education. Over the past year, this project has been 
successful in attracting many students over and above those attracted by financial 
support. The u* students, who participate on the project in this fashion, not only 
gain academic credit but also share in the challenge and the stimulation offered 
by the program. This past year is no exception with a total jt 37 students actively 
involved in the project, only 12 of whom received some financial support. Of this 
total, two were doctoral level, twelve were master's level and twenty-three were 
undergraduate level students. 


The students' reactions regarding the educational value continue to be 
very positive as evidenced by the number: of new students attracted each year. The 
f culty appraisal of value is also very high. The project provides a unique learn- 
ing experience in an academic setting for a number of reasons. These and other 
aspects of the educational impact arc discussed in detail in Appendix A. 

First, all of the problems addressed by the students are relevant and 
significant. The student realizes that the problem must be solved within the con- 
straints that apply anc not within artificial constraints which would make the 
solution relatively straightforward. Faculty attitudes and guidance arc very im- 
portant in this regard as are those of the technical monitors from Jet Propulsion 
Laboratory whose role is described below. 

Second, the project is sc organized and administered as to create a mean- 
ingful Industrial clinrte. The students are aware of the short and long-term techni- 
cal objectives nnd of their role in tne project. The problem of interfacing their 
task with those tasks on which it impinges becomes painfully clesr. fhey realize 
that interfacing with other tasks is just as crucial a constraint as are the direct 
technical objectives. The students are expected to plan their task, to formulate a 
milestone schedule and to achieve their stated objectives. Even though the students 
tend to be overly optimistic, they lean: from this experience how to go about plan- 
ning and executing their assignments. The studenls are required to make periodic 
progress reports and to prepare such progress and technical reports as arc n ccssary. 
All in all, the students are experiencing the responsibilities of practicing engineers. 
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Third, progress review meetings involving the project technical monitors 
from JPL and invited guests are held at approximately ten week intervals. These re- 
view meetings are one of the most vital aspects of the project both from the techni- 
cal and educational points of view. Students make presentations of progress to the 
group of monitors, guests, faculty and the other participating students. The progress 
is reviewed critically by all participants with appropriate comments and discussion. 
These meetings develop the student* s ability to prepare his materials, including 
visual aids and demonstrations, to deliver his presentation and to react to questions 
and suggestions. The monitors and guests frequently provide additional information, 
knowledge and guidance which promotes the progress of the work. The progress review 
meetings are of considerable value in maintaining steady and coordinated progress. 

* | 

In summary, the project has contributed substantially to the education of 
engineering students and the professional growth of faculty interested in advanced 
study in the area of desif u 
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V. Detailed Summaries of Progress 

Task A. Vehicle Conf iquration , Control, Dynamics, Systems and Propulsion 

The objectives of this task are two- fold: (a) to design, construct and 

evaluate a data and sample gathering rover to augment a second generation Viking 
lander or to enable a sample return mission for an extensive Mars exploration, and 
(b) to provide a tost, bed for the evaluation of alternative hazard detection and 
path selection systems for an unmanned autonomous rover. 

Towards the accomplishment of these objectives, this group has addressed 
the following areas: vehicle dynamics, vehicle structure and propulsion, wheel 

design, and telemetry and control systems. These efforts are described i: more 
detail below. 
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Task A.l. Experimental and Analytical Stu d ies of Vehicle Dynamics - P. Marino 

Faculty Advisor: Prof. G. N. Sandor 

The objective was to analyze the vehicle's dynamic behavior and modify 
the vehicle's suspension so as to improve stability and reduce shock on the payload. 

A two degree of freedom model of the vehicle's dynamics was developed and 
was used to choose appropriate rotational dampers for the rear wheel struts, Figures 
2 and 3. These dampers have non-linear characteristics, and a mathematical model 
of the dampers was incorporated into trie dynamic analysis program. It was found that 
the dampers must be geared to the rear torsion bar in such a way as to reduce the 
torque load on the dampers and increase their rotational velocity. A chain drive 
reduction was chosen because of light weight, high torque capability, and relative 
immunity to dust and foreign particles. A ratio of 6 to 1 was chosen and an appro- 
priate chain drive installed, li addition to implementation of rotary dampers, it 
was found that a slightly stiffer torsion bar would improve the vehicle's dynamic 
behavior. A suitable torsion bar was designed, fabricated and installed on the 
vehicle model. Future work includes dynamic testing of the improved vehicle, model so 
as to verify the mathematical model and analysis. 


Task A. 2. Vehicle Structure - J. Harrison and B. Musits 
Faculty Advisor: Prof. C. N.* Sandor 

This task involved: redesign of the front end of the vehicle to overcome 
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structural defects and weakness observed during the field testing; implementation 
of front wheel drive; implementation of strut sensors to provide vehicle configura- 
tion data; latch sensors to insure locking of members on completion of deployment; 
redesign of the torsion bar system and payload box; and construction of a laboratory 
obstacle course. 

Front End Structural Revisions 

Field testing of the vehicle and subsequent examination of 
components indicated that a number of structural improvements 
to strengthen the components were required. These have been 
made and results to date are satisfactory. The ability of the 
front end to withstand shock has been increased substantially' 


Front Wheel Drive 


Field testing indicated that the anticipated obstacle negotia- 
ting capability of a vehicle relying solely on powered rear 
wheels did not materialize. Analytical studies indicated that 
a significant improvement in step climbing ability could be ex- 
pected with the addition of front wheel drive, Figure 4. A 
front wheel drive system was designed, constructed and implemented. 
Preliminary testing indicates that front wheel drive will add sub- 
stantially to the vehicle’s ability to negotiate potentially 
hazardous terrain situations. The vehicle’s ability to overcome 
step obstacles depends not only on the torque applied to the wheels 
but also the traction between the wheel and the terrain. The new 
wheels described under Task A. 3 are expected to provide a sub- 
stantial increase in performance beyond that shown in Figure 4. 

The addition of front wheel drive has also resulted in an improved 
reliability of the deployment maneuver. 


Strut Sensors 


Successful autonomous roving requires not only an effective 
hazard detection system but also knowledge of the state of the 
vehicle such as pitch and roll. For example, a terrain feature 
which may be entirely negotiable if the vehicle is essentially on 
a horizontal plane may be a severe hazard if the vehicle is climb- 
ing a 15° slope. In order to provide the necessary vehicle ’’state" 
data, sensors which can detect the orientation of the rear and 
front struts relative to the payload have been installed along with 
sensors detecting the inclination of the front axle as well as the 
steering angle. Arrangements have been made, Task A. 5, to trans- 
mit the sensor data to the remote control station. 


Latch Sensors 

The launch configuration requires the unfolding, extension and 
latching of both the rear and front struts. Four latch sensors, 
one for each rear wheel strut and two for the front strut were 
installed. Arrangements have been made to transmit the signal that 
the latches have been driven home. 
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WITH FRONT WHEEL DRIVE 
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SLOPE (DEGREES) 


Fig. 4. Effects of Adding Front Wheel Drive 
on e Combined Slope-Step Obstacle 





Torsion Bar and Payload Box Design 


The torsion bar and gearbox were redesigned to eliminate 
observed structural weaknesses and to permit modification of 
the payload box- The previous design. Figure 5, can be com- 
pared with the current design, Figure 6. Not only has the 
torsion bar gearbox and assembly been improved structually but 
its reduction in size has increased the payload volume available 
for necessary on-board systems and reduced the payload weight by 


40 %. 


Laboratory Obstacle Course 


Although field testing is the ultimate test site which can re- 
flect the realities to be overcome by the vehicle, it does not 
provide the quantitative data required for an in-depth appraisal 
of the vehicle's capabilities. Such an appraisal can lead to 
design changes which can result in significant improvements ir 
vehicle performance. To this end, a laboratory obstacle course 
has been constructed. This system, Figure 7, is very flexible 
and provides the project group with the ability to create special 
terrain features according to the indicated needs. The present 
course which measures some eight feet by sixteen feet can be ex- 
tended without difficulty. 


Task A. 3. Wheel and Grouser Analysis and Design - P. Marino 

A. 3. a. Wheel Tester and Grouser Design 

The objective was to find an optimum tread for use on the RPI MRV 
vehicle model while examining grouser behavior in general. 


Last year's tread design, while yielding good traction, proved to lack 
durability. A more durable tread material, commercially available for 
use on conveyor belts, was tested and chosen for use on the vehicle model. 
This tread material was applied to the stiff-rim type wheels in use last 
year and yielded satisfactory results. Current changes in wheel design 
at this time required a wrap-around type tread with reasonable backing 
stiffness to be used on a soft-rim design wheel. A search was then con- 
ducted for suitable tread material. The best readily available grouser 
was found to be the outside rim of an all-terrain-vehicle type motorcycle 
tire. One of these tires was obtained and modified for use on the new 
soft-rim type wheel. Future work includes testing of this new type grouser 
on the RPI wheel tester to determine its effectiveness on varied surfaces. 



A. 3.b. 


Wheel Analysis 


The objective was to analyze the existing RPI MRV model wheels and decrease 
footprint pressure while maintaining lateral stability and Increasing trac- 
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In an effort to reduce footprint pressure a change to a soft-rim design 
was investigated, Figure 8. It was found that a split hoop design combined 
with a soft rim yields adequate lateral stability, low footprint pressure, 
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OBJECTIVES: TO MAXIMIZE TRACTIOti KID TERRAIN MOBILITY 

DECREASE FOOTPRINT PRESSURE WHILE RETAINING 
j LATERAL STABILITY 



Fig. 8. Old and New Wheel Design 
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and high load capacity before overloading of the wheel hoops. This 
new design incorporates an inner split hoop and an outer grouser material. 
The soft outer hoop yields under low load and causes a large tread area in 
contact with the ground and thereby a low footprint pressure. Under higher 
loads, the outer hoop contacts the stiffer inner supporting steel and the 
'the stiffer springs carry the additional load. In this way, the wheel 
combines two desirable characteristics: a low footprint pressure, and the 

ability to handle high loads due to the inner hoops. Figure 9. 

The new wheel design offers a distinct advantage in the negotiation of 
hazardous terrain in that its softness will provide a better grippage and 
increased traction over the past steel rim wheel, (see Task A. 2 Four 
Wheel Drive). It is hoped that step climbing ability will far exceed 
predictions for the stiff wheel. 

Future work includes optimization of the design in terms of hoop stiffness, 
number of hoops, and spring hoop arrangement. 


A.3.C. Detail Wheel Design and Construction 

The objective was to combine the wheel structure and available grouser 
material and construct wheels for the RPI MRV vehicle model. 

The new double hoop wheel design has been combined with the tire-type 
grousers and a new type wheel has been constructed. Plans include test- 
ing of this completed wheel and construction of four such wheels for use 
on the RII MRV vehicle model. 


Task A. 4. Telemetry Systems and Vehicle Electronics - T. Barren, J. Cooley, 

K. Fell, T. Geis, T. Kasura, R, Richter, R. Pfeiffer 
Faculty Advisors: Profs. G. N. Sandor and D. Cisser 

Attention this year was focussed in three major areas- (a) update and 
modifications of the vehicle electronics, (b) implementation of front wheel drive, 
power and control systems, and (c) the construction of a data link trom the vehicle 
to an off-board computer. 

Modifications to Existing System 

A new control system, found to be necessary during the testing of the 
previous summer, was designed and built this year. This control box. 

Figure 1C, facilitates more rapid commands to the vehicle and a more 
effective human-machine interface than the previous design and hence more 
precise control of the vehicle. The control box is adaptable to inter- 
face with an offboard minicomputer to allow control of the entire system 
when under closed loop control. 

Modifications were also made to the onboard vehicle electronics. A 
problem encountered during previous vehicle testing consisting of an 
overload in the steering drive circuitry was corrected by changing 
the function of steering motor control relays. Voltage fluctuations 
in the vehicle electronics power supply were found to be the cause of 
"creep" In the rear motor drives and variations in motor speed. To 


! 



Fig. 9. Deflection Characteristics of the Double Hoop 
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minimize these fluctuations, a commercial power supply, capable of 
supplying 1 amp at both plus and minus fifteen volts with an input of 
twelve volts, was used instead of the unit built previously, which had 
been operating a*, near maximum ratings. This new power supply is work- 
ing at about 1/3 capacity, allowing for a great deal of expansion while 
still maintaining well regulated output voltages. 

Adoption of the mini rover concept this year necessitated the capability 
of much lower vehicle speeds in order to acconmodate the laser hazard 
detection system, Task B, as opposed to the laser rangefinder system, 

Task D.l. To obtain these slower speeds, the motor drive circuitry was 
modified by adding gain to the feedback stage to allow the drive motors 
to turn more slowly. 

Other modifications included a sensor to inform the vehicle electronics 
when the vehicle was in the "flipped" condition. This data was necessary 
to invert the differential action on the rear wheels to allow the vehicle 
to be driven in this mode. Mast control circuitr vas added to allow a 
roast to be raised and lowered by remote control. 


Front Wheel Drive 


Electronic control for front wheel drive was implemented this year as a 
simple parallel or series connection over the main drive batteries. The 
parallel connection provided full power to the front wheels while the 
scries connection provided a differential action to allow steering over 
large angles. This system provided adequate mobility but also showed 
that a four wheel speed control system would he desirable in terms of 
power consumption, precise speed control and steering at a standstill. 
Under consideration at this time is a four wheel speed control system 
that will achieve these benefits as well as have a fourfold increase in 
the accuracy of each wheel's speed as a function of steering angle. This 
new system is scheduled to be implemented by the middle of this summer. 

Vehicle to Offboard Computer Data Link 

The mai:. objective this year was to send information concerning the 
state if the vehicle to an offboard computer. This is intended au a 
partial implementation of Phase II shown in Figure 11. Phase III which 
represents a fully autonomous, closed loop control of the rover is 
scheduled for implementation by March 1976. 

Located on the vehicle are four latch sensors to verify the deployment 
of the vehicle, four tachometers to monitor the speed of each wheel, seven 
potentiometers to measure the positions of the front axle, emergency 
maneuver motor, steering angle, elevation of the front and reai strutc, 
deflections in the torsion bar, and a gyroscope to monitor the pitch and 
roll of the vehicle. Figure 12. 

Each sensor is sampled, cc /ertc-i Into a 12 bit digi:al word, and pre- 
sented with the latch sensor data to a digital multiplexor which selects 
this vehicle data or the data from any hazard detection system mounted on 
the vehicle. The multiplexor is capable of taking this vehicle data or 
hazard detection data any any desired "duty cycle" or having a "hazard 
detector override". 
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Fig. 11. Steps to s Completely Autonomous Ro ■* . Vehicle 






Fig. 12. Vehicle Attitude Sensors 


25 . 


\ 


A digital "word" is constructed which consists of a start bit, an 
address which locates the data word's origin to the hazard detector 
or the sensors, and if from the sensors, which one ( The data is then 
put in serial form, followed by an error correcting code, four st^p 
bits and one framing bit, Figure 13. The error correcting code is a 
combination of a five bit Hamming code and one parity bit. This code 
makes it possible to correct any single error and detect any multiple 
error. The system is currently sending 500 words per second and is 
capable of sending much more, if necessary. 

The receiver unit, Figure 14, located at the computer, recognizes the 
start bit as the first zero received, and starts to count off the rest 
of the word changing it from serial to parallel. The start bit is 
verified to insure that the receiver didn't trigger on a noise pulse 
and finally, the four stop bits are verified to insure the first zero 
was indeed the start bit. If all is correct, the word is stripped of 
start and stop bits, any single errors are corrected (words with multiple 
errors are aborted), the word is fed into a register and a "data ready" 
pulse is created for the computer to use. 

The lowpower telemetry transmitter uses direct carrier frequency shift 
keying in the VHF band, Figure 15. It is physically small and designed 
to handle at least 50 kilobits per second. It consists of two stages, 
a high frequency voltage controlled oscillator, fed by the serial digital 
data which drives a lower power output stage. 

The receiver is a commercial unit chosen for its high overall bandwidth, 
Figure 16. It consists of a radio frequency amplifier, oscillator-mixer 
stage, and intermediate frequency amplifier section with an overall band- 
width of four megahertz. An FSK, Travis type detector was designed to 
derive the data from the received signal. Variations in carrier frequency 
are interpreted as digital symbols allowing the detector to be quite in- 
sensitive to amplitude variations and noise. 

A display was built to show the data from the vehicle that samples the 
data at the computer interface. It displays at all times the state of 
the latches on light emitting diodes. In addition, the operator can 
select any sensor on the vehicle and watch the transmitted information 
from that sensor on a digital readout. 

Also this year, preliminary work was done on a computer program which 
will use the information from the vehicle to create a graphical display 
of the state of the vehicle. This would show top, front, and side views, 
as well as the roll and pitch on the CRT display, Figure 17. This will 
assist operation of chc vehicle from the computer. 


Task B. Hazard Detection Systems - S. Cairns and D. Holly 
Faculty Advisor; Frof. D. Gisser 

The Mars Roving Vehicle must be capable of operating autonomously over 
varied terrain for extended periods, although a reasonable number of (automatic) 
stops with requests for human assistance are to be expected. Optimum operation 
would include continuous scanning of the terrain in the general direction of motion 
in order to select a "best" path and initiate avoidance maneuvers early. Systems 
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Fig. 13. Transmitter Encode Logic 


















Fig. 14. Receiver Decode Logic 



















Fig. 15. RPI MRV Telemetry Transmitter 







Fig. 16. RPI MRV Telemetry Receiver 





(Arro - indicate changes in vehicle position revealed in that particular view) 
Fig. 17. Typical IDIIOM CRT Display 




of this sort were planned during the past several years, and they require a range 
finding sensor with a range resolution of a few centimers ( 5 cm) at ranges to 
30 meters, appropriate angular resolution, and scanning capability. Such range- 
finders can be built using pulsed laser sources, multiplier phototube detectors, 
and optical systems with an aperture of 5 to 10 cm. Implementation of such rangers 
with solid state detextors is just about possible theoretically and should become 
feasible in practice with advances in both laser and detector technology. The wide 
aperture required, and the necessity for a mechanical scanning arrangement makes 
the outlook for a system with the relatively short scan time required, somewhat 
dubious. Such a system, however, is still considered a long term possibility 
(Task D) and improvements in electro-optical materials may make a solid state scan- 
ning system feasible. 

The short term objective has changed. Instead of a precise + 5 cm, long 
range 3-30 meters, time-of-f light system, whose basic purpose is to supply detailed 
terrain data in quantity, the current objective is to impelement a short range, 

(0.5 - 2.5 meters) low resolution, triangulation system whose basic purpose is to 
supply the vehicle with information of the go-no-go variety. As a minimum, it is 
expected that both boulders and holes larger than some preset threshold can be 
detected with this arrangement. It is also possible that laser scanning and multiple 
detectors may provide rather detailed information regarding the terrain in the 
immediate vicinity of the vehicle. 


Principles of Operation 

The system is basically a reflected light detection arrangement, 
consisting of a light source with a relatively narrow beam angle 
directed at the terrain at a point of interest about one meter 
in front of a wheel. A light detector with a somewhat larger 
beam angle is located at a point as far as is practicable from 
the source, as shown in Figure 18. Reflections from the ground, 
with or without small variations in level, are received continu- 
ally and are considered nonnal. Any very sudden drop off or 
barrier, however, will result in no reflections within the 
acceptance angle of the detector, as shown in Figure 18. This 
system for the detection of obstacles and drop offs is not new, 
but in actual use for a similar purpose as a mobility aid for 
the blind. 


Implementation 

A continuous source of light with high enough efficiency for 
high output with low power requirements is not available, but 
the injection laser can be pulsed to a very high output for a 
short time, making its average power drain quite low. It has 
the additional advantage of producing its light output from a 
very small area and over a small bandwidth, permitting a simple 
optical system to produce a narrew beam. The emission is around 
9000 A, in the near infra red, a wavelength where semi-conductor 
detectors are quite sensitive and most minerals have reasonably 
nigh reflectivities. The transmitter circuit is shown in 
Figure 19. Operating from a supply source of 15 volts, it uses 
a total of 150 milliwatts to produce about five hundred 150 
nano-second light pulses per second. It is packaged in a box 
approximately 14 r 10 x 9 cm, supported by an adjustable ball 
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Fig. 18. Hazard Detection System Concept 








and socket mount. The lens is 8.75 mm in diamter and has a 
focal length of 35 mm. Notice that there is an internal dc 
to dc converter to charge the storage capacitors to about 80 
volts. 

The receiver circuit is shown in Figure 20. It consists of 
a 25 m diameter lens with 21 mm focal length focusing light 
on a 5D PIN photodiode, followed by a two stage low noise 
wide band preamplifier operating from a nine volt battery. 

The total power drain is 32.5 mw. The receiver box is about 
10 x 7 x 5 cm. which is also ball and socket mounted for easy 
adjustment. The output feeds an amplifier and pulse stretcher, 
shown in Figure 21. Its output will be eventually rectified, 
and the absence of a signal used to indicate the presence of 
an obstacle or hole. For experimental work, however, it is 
very convenient to have an audible indication of reflections 
so the circuit shown in Figure 22 is used. Strong reflections 
produce a continuous audio tone. Marginally weak reflections 
produce an intermittent, raspy tone, and no reflections, no 
tone. 


Discussion 


The system has been shown to detect obstacles easily and drop- 
offs of about one third meter at about a two meter range with 
a one meter baseline, although echoes are weak from some very 
dark materials. Work is progressing on more carefully defin- 
ing the receiver beamwidth needed and in considering the 
relative advantages of several configurations of several de- 
tectors and sources in order to better protect thr vehicle. 


Task C. Path Selection System Simulation, Development and Evaluation - 
D. Matthews, D. Sharp, D. Campbell, J. Krajeweski 
Faculty Advisor: Prof. D. K. Frederick 

An autonomous Mars Roving Vehicle must detect and avoid obstacles on its 
way to some predetermined target. This task is aimed at the development and evalua- 
tion of sensing and avoidance schemes using a digital computer simulation. In the 
past year, progress has been made in the development and evaluation of both mid- 
range and short-range path selection systems. Improvements have also been made on 
the computer simulation package; improvements which will make the package a more 
effective design tool for the new project goals. 

A mid-range path selection system with two and three beams fixed on the 
vehicle has been developed and tested. Uniform and non-uniform axlmuth scans have 
been employed and encouraging results with positive obstacles have been obtained. 

In addition, a short-range system has evolved from a basic terrain modeling system 
with simple go/no-go conditions into a more sophisticated path evaluation scheme 
with memory capability. Further, the short range system has demonstrated its ability 
to avoid craters. Improvements have been made in the group's digital simulation. 
Including modifications in the display of output so that results are more easily, 
read and interpreted. The terrain generation and vehicle dynamics simulation are 
being improved so that tho evaluation of the newly proposed triangulation hazard 
detection system will be easier and more realistic. Work is also underway to 







Fig. 21 . Amplifier and Audio Converter Electronics 





simulate the new sensor so that the triangulation scheme can be included in the 
digital simulation and an appropriate path evaluation algorithm can be developed 
for the vehicle. 


C.3.a. Development of a Mid-Range Path Selection System 

The objective of this subtask is to investigate the feasibility of a 
vehicle-fixed terrain scanning concept for mid-range path selection 
purposes. An optimal sensor configuration and path selection method may 
then be determined for this type of sensing system. 

Analysis of a single beam sensing system, Ref. 1, has suggested refine- 
ments which have been incorporated in succeeding multi-beam systems. 

Both double and triple beam systems have been proposed and evaluated. 

The most recently evaluated system is an advanced, three beam system 
which employs nonuniformly spaced azimuth scans. 

The general multi-beam system is described below and the results of 
the evaluation of the double beam version are briefly discussed. The 
advanced three beam system is then described and discussed in detail. 

Multi-beam System Description 

In general, a path selection system must sense the terrain, model the 
terrain using data from the sensor, and using the terrain model, select 
a path. Described below are the three distinct subsystems that perform 
these functions: the sensor, the terrain modeller, and the path selec- 

tion algorithm. 

Sensor Systems - The computer will simulate up to four sensors, 
mounted on a vertical mast which is fixed to the front of the 
vehicle. One can specify the number of sensors to be used and 
tneir heights on the mast. Each sensor scans in a plane per- 
pendicular to the mast, making twenty-nine range measurements. 

The azimuths at which measurements are made are also specified 
by the user and can be uniformly or nonuniformly spaced. Note 
that this sensing scheme works for positive obstacles only 
(i.e. those obstacles which protrude above the ground). For a 
negative obstacle, an ideal sensor is simulated which senses 
the obstacle after the vehicle bumps into it. 

Terrain Modeller - Taking the twenty-nine range measurements 
from the sensor, the terrain modeller generates a go/no-go 
map for fifteen possible forward paths. This is done by com- 
paring the range measurements with a computed set of thresholds, 
which take into account the vehicle width and a specified 
buffer zone. 

Path Selection Algorithm - There are two modes of operation 
for the algorithm. In normal operation, the algorithm sets 
the vehicle on the most direct, acceptable path to target. 

When all possible fifteen forward paths are blocked, the 
emergency mode is implemented. First, the vehicle is backed — - 





up in a straight line. Then a new sensor scan is made and 
the most direct, acceptable path to target is again taken, 
excluding the seven forward paths which would again lead the 
vehicle across unsafe terrain. The emergency mode can also 
employ a memory capability at the user's option. If memory 
is used, a preferred side is chosen after taking a new sensor 
scan. This choice is made by summing the ranges on each side 
of the vehicle and selecting the side whose sum is largest. 
For a specified number of subsequent scans, the terrain 
modeller blocks paths on the unpreferred side. 


Double Beam System Evaluation 


The multi-beam mid-range system was first investigated using two sensors 
mounted on the mast at heights of 0.5 meters and 1.5 meters above the 
ground. This system was evaluated using the standard testing procedures 
described in Ref. 2, namely. 


Case I 2x2 meter boulder encounters 

on flat base terrain 


Case 

II 

2/3 x 2/3 meter boulder 
on flat base terrain 

encounters 

Case 

III 

Gently rolling terrain 

/ 

Case 

IV 

Rolling incline 


Case 

V 

Boulder-crater field of 
complexity 

moderate 


Each case was simulated in the absence of noise and witii the addit - ;f 
attitude and range noise. Combinations of range and attitude no , ” * 

also simulated. 

During the simulation, the obstacle avoidance bull rs w kept t 0.5 

meters for the lower sensor and 1.0 meters for t i,e u pro, senior 4 nal’ r 

buffer on the lower sensor resulted in a lower 'e. ^ieivity t( . # itude nc \e 

and prevented the vehicle from overreacting to boulders. y . iv.ng 

the higher sensor a larger buffer, the vehicle had ' >. ,h tim p - ; *i ' : L 
avoidance maneuvers around larger obstacles, avo: 1 r w mnece^ary back ps. 

As indicated in Fig. 23, the system was capab of avoiding obstacles of 

a size comparable to the size of the vehicle tut as the size 

boulder decreased, the system hd \ difficu. \ v ^ing th* i ± ^ for avoid- 
ance maneuvers to be conveniently employe . i :ly, the system lost 

sight of the obstacle between successive beam shots indicating a need for 
decreased azimuth differences. This : ^crea °.d r imuth difference was in- 
cluded in the three beam system, vLlc.i is uescti.bed later 

As a result of wandering noted or> rot Ling and inclined terrain, the thro*. ->ld 
on the lower sensoirwas reduced, a a nonui.iformly spaced differ* s 
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was used. The configuration which was implemented is shown in Fig. 24. 
Since the azimuth differences decrease towards the front of the vehicle, 
the vehicle was not able to turn as sharply as before. Eut with de- 
creased azimuth differences at the front of the vehicle, it was thought 
that the vehicle might not have to turn sharply since the vehicle could 
start avoidance maneuvers well in advance of any obstacles. Two possible 
alternative approaches were proposed for better vehicle performance on 
turns. One possible approach would be to slow the vehicle down when a 
sharp turn is required. Alternatively, one could increase the number of 
beam shots per scan and the number of possible paths. To avoid undue 
complication in the path selection algorithm, the former proposal was 
implemented. 


Three-beam Mid-range System 

As recommended by the two-beam study, another sensor was added at 1.0 
meter, allowing for a lower threshold on the lower sensor, Ref. 3 and 4, 
Also, the nonuniform scanning scheme was carried over from the two beam 
system to increase the information known about the terrain in front of 
the vehicle. As indicated in Fig. 24, azimuth spacings were decreased 
towards the front of the vehicle. 

It was noted that the three beam sensor system, with decreased thresholds 
on the lowest sensor and nonuniform scanning, wandered less than its simple 
two beam predecessor, as evidenced by comparing Fig. 25 with Fig. 26. It 
was found that on hilly and inclined terrain, the target was reached with 
a minimum of wandering about, but that there were a number of unnecessary 
backups required for small boulders on flat terrain. There appears to be 
a trade-off between performance with single boulder encounters and perform- 
ance on an incline. Decreasing the thresholds for reduced wandering on 
inclines led to unnecessary backups with single boulder encounters. It is 
recommended that the thresholds be determined by vehicle attitude. In- 
creasing the thresholds with a smaller in-path slope could lead to fewer 
backups on small boulders in flat terrain. As vehicle attitude increases, 
thresholds can be decreased so that the vehicle will not wander and will 
head directly for the target. 

Simulation results on complicated bouluer-crater fields, Fig. 27 and 
Fig, 28, show the trade-off that was made, Note that improved performance 
on the sloping terrain has resulted in numerous backups and a path whose 
figure of iueric is low due to the lower thresholds necessary to prevent 
wandering on an incline. Further, the low figure of merit resulted from 
the mid-range system's inability to detect negative obstacles. 


C,3.b. Short Range System 

The objective of this subtask is to simulate a path selection system 
having sensing capabilities limited to a range of three meters, Ret. 5. 
It is hoped that the simulations run in this investigation will be 
sufficient to decide if the short range system can be used as: 

1) a primary guidance system, 

/ 

2) a guidar ce subsystem used with a midrange 
subsyst * to detect small positive obstacles and 
craters, and/or 

3) a backup guidance system. 
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fiG. 26. Triple Bean System vith Reduced Thresholds 

on the Lover Sensor and Hon-Uniform Scanning 
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Note that l*h * censor described below is a conceptual configuration 
only and i* rut intended to indicate the need ior a specific hardware 
design. 

As in the mid-rorge case, the path selection system consists of a 
sensor, a terrain modeller, and a path selection algorithm. Each is 
described below. 

Sensor - The sensor, analogous to a blind man with a cane, 
consists of two proximity detectors mounted at the end of 
an arm -attached to the vehicle. A picture of the system 
is shown in Fig. 29. The arm sweeps from side to side and 
takes measurements J n the forward and downward directions. 

Forward measurements were used to protect the arm by retract- 
ing it. Under normal operating conditions, the arm is two 
meters long, but it can bo retracted to one meter if necessary. 

The downbeam measurements provide readings of the height of the 
terrain at the end of the arm relative to the height of the 
vehicle. This information is used in the terrain modeller. 

Terrain Modeller * The terrain modeller makes use of three 
successive sweeps of the vehicle's arm to fit the terrain 
in > ach of the possible vehicle paths to a height function 
which i vtidratic in the horizontal and vertical distances 
from the vt-iiicle. Each path, illustrated in Fig, 30, is then 
evaluated by a path selection algorithm. 

•_v jc- : ;n 1 gori t hm - The short range system takes the 

-.rect : argot which meets all of the following criteria: 

1/ the in-path slope does not exceed the vehicle's 
ibility to traverse it. 

2) he i -path slope does not exceed the maximum 
ogle at which slippage is bound to occur. 

N the change in both cross-path anil in-path slope 
is such that the vehicle will not be hung up or 
"ill hol "bottom out ’. 

4) there art no positive obstacles in the way, as 
found from the forward proximity sensor. 

oim i ation Result s 

p relimxiM*ry testing was run on the algorithm, following the standard 
testing procedure that was previously mentioned. Single boulder and single 
crater runs were made, t> encouraging results. In Fig. 31, we see that 
the proximity s isor s _>fully avoided a positive obstacle using the 
i^rv~rd roxiraity n aid the arm's capability to retract. Although 
there as some vanui h.c, obstacle avoidance proceeded nicely and the target 
was reached 


/ 


a) Top View 



b) Side View 



Fig. 29 - Short-Range Systcu Sensor Configuration 
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Selection in the Short-Range System 
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A single crater encounter was simulated which resulted in this group's 
first success ul negative obstacle avoidance maneuver- The simulation 
is shown in Fig. 32. 

One further feature of the system's path selection algorithm is the 
memory capability. Conceptually, information from the successive scans 
is stored so that when all paths are blocked, the vehicle can return to 
a previous decision point and find a new path to target. This concept 
is illustrated in Fig. 33. On scan 3, the vehicle decided to take a 
path around the boulder but later on scan 6 found itself in impassable 
terrain. The vehicle can then return to the decision point and circum- 
vent the boulder by going around the opposite side. An actual simulation 
of this capability occurs in an early run of the algorithm, in Fig. 34. 
The vehicle found itself proceeding into the crater, having erroneously 
selected a path towards the crater from a previous scan. At the position 
marked X, the vehicle models the terrain and discovers that this path 
would not be passable. Finally, the vehicle uses its memory capability, 
returns to the decision point and moves around to the target. While the 
capabilities of the algorithm have not been fully tested, these initial 
insults demonstrate that the short range concept warrants further study, 
refinement, and simplification. 


General Concl u sions 

In viewing the work done to date, certain conclusions of a general nature 
can be made which should be applicable to sensing schemes other than the 
particular ones studied. Work done on the mid-range path seleccion system 
has led to studies of sensor placement, threshold distances, nonuniform 
scanning, and the effect of vehicle speed on the vehicle’s turning ability. 
Vehicle memory has proven to be helpful when the vehicle is confronted with 
an emergency situation. Some general conclusions may also be drawn con- 
cerning the terrain modelling schemes that the Path Selection Systems Group 
has employed to date. 

In the mid-range system (a system which detects positive obstacles on the 
order of ten meters away), one, two and three beam systems have been studied. 
On multiple beam sensing systems, the lower beams employ a lower threshold 
than the upper beams since small rolls and pitches can cause the lower beam 
to sense an obstacle when there are none. Also, the upper beams require 
longer thresholds so that obstacle avoidance of larger obstacles can begin 
with enough time for the obstacle to be successfully avoided without neces- 
sitating emergency backup maneuvers. 

While the two beam system was being tested, it was found to have difficulty 
avoiding boulders which were smaller than the vehicle, but big enough to 
block vehicle passage. These boulders tended to be missed by the sensors 
because it was between the aximuth scans. Hence the scans were concentrated 
towards a more limited field of view, but most of the obstacles with which 
the vehicle must concern itself tend to be in front of the vehicle. Im- 
proved path selection performance resulted from these nonuniformly spaced 
azimuth scans. 

After centering the scans twoard the front of the vehicle, the vehicle was 
found to be unable to turn as sharply as was previously the case. The 
reason for this was that the paths to be selected were also centered towards 
the front of the vehicle. Two possible solutions were proposed for increas- 
ing the vehicle's ability to turn. First, the beam shots per scan could be 
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Fig. 33. Effects of Vehicle Memory for the Short-Range 
System 
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increased, thus increasing the number of paths available to the vehicle. 
Alternatively, vehicle speed could be reduced in the turn, making a 
wider range of paths ‘available to the vehicle. The latter proposal was 
implemented since it did not complicate the path selection algorithm. 

When ail paths to the vehicle are blocked, the emergency mode of the 
path selection algorithm is activated. The vehicle must remember where 
the obstacle is or it may tend to bump into the obstacle repeatedly. 

Hangups of this sort were avoided using ’ chicle memory to block off paths 
which would again lead the vehicle towards the obstacle. This kind of 
algorithm is necessary in virtually any emergency mode operation of the 
path selection algorithm. Another possibility, especially with a short- 
range type of system, is to return to a previous decision point, a point 
where one path was chosen over another due to better target proximity. 

The alternative path can then be explored, decreasing the chances that 
the vehicle will be boxed in. 

In modelling the terrain in front of the vehicle, computing the in-path 
slope from the height of the vehicle and the height of the short range 
sensor leads to an erroneous model of the actual terrain. That is, it 
assumes that the terrain is uniformly sloped between the sensor at a ' 
known arm length away, and the vehicle. Algorithms which employed this 
type of modelling performed poorly. An approach which led to the detec- 
tion of a negative obstacle was to model the terrain using three previous 
sensor scans. A quadratic fitting algorithm was then used, and the cross- 
path and in-path slopes were then approximated. Although this algorithm 
has not been completely tested, preliminary testing indicates that negative 
obstacle avoidance is possible. 

In light of the slower speeds and reduced terrain size of the new mini- 
rover concept, more realistic terrain representations and vehicle dynamics 
simulations have been formulated. These improvements are being incorporated 
in the simulation package. 

The triangulation sensor proposed under Task B is being modelled and will 
be simulated so that appropriate path selection systems based upon this 
concept can be simulated and design alternatives studied. This work will 
provide direct support for the ongoing vehicle and sensor hardware work 
during the coming academic year. 

Documentation of the simulation program will continue to receive a con- 
siderable amount of effort in order to improve its effectiveness. 


Task D. Obstacle Detection and Terrain Modeling 

Investigations continued in the areas of obstacle detection by edge per- 
ception, computer graphic' evaluation of edge detection procedures, terrain modeling, 
range measurement scanning concepts, laser range finder errors and obstacle 
classification. Summary progress reports covering each of these topics follow 
below. 



Task D.l.a. Obstacle Detection - Ranjan Sonalkar ' 

Faculty Advisor: Prof. C. N. Shen y 

\ 

The objective of this task is to develop a practical and reliable obstacle 
detection scheme which will make use of the range matrix obtained by a 
laser rangefinder. The range measurements may contain errors of anywhere 
between + 5 cm to + 10 cm. The Four Directional Ratio (FDR) algorithm 
developed last year has been documented in Ref, 6, The algorithm takes 
advantage of the fact that a significant change in the range readings occurs 
as the laser crosses the obstacle edges. For example, as the laser crosses 
the top edge of a boulder (positive obstacle) or the near edge of a crater 
(negative obstacle), the measured range increases suddenly and the FDR; 
identifies these edges by detecting the corresponding range changes. How- 
ever, th° bottom edge of a boulder and the far edge of a crater, constitute 
changes in slopes and not ranges. Ihe change in range caused by the change 
in slope, is then too small to be detected by the FDR algorithm. However, 
complete characterization of the obstacle is necessary for the rover to 
successfully negotiate the path towards its destination. Therefore, it was 
deemed necessary that the FDR algorithm be improved or a new one designed, 
which will fulfill the requirements. 

The new algorithm that is proposed and is an adaptation of the Rapid Estima- 
tion Scheme (RES) described in Refs. 9 and 10 is documented in Refs. 7 and 8. 


Outline of the Scheme 


The range measurements contain errors, presumably of known statistics. There- 
fore, the measurements are not exact and cannot be fully relied upon. Also, 
the slope changes at the edges of obstacles have small effects on the range 
measurements and the effects are of the order of magnitude of the errors. 
These difficulties suggest the use of some kind of a reliable probabilistic 
mechanism to detect the edges, which will also be able to distinuish between 
the measurement errors and the changes of slopes at obstacle edges. 

The scheme processes each column of the range matrix to determine whether 
an edge intersects it, Fig. 35. Each column is processed upwards, element 
by element and three consecutive elements are considered simultaneously. 

Three hypotheses are defined as follows: 

Edge exists between stage i and i+1 
2 : Edge exists between stage i+1 and i+2 

Edge, if any, exists after stage i+2 


The three hypotheses are represented by the 3 branches of a decision tree, 

Fig. 36. Decision theory is employed to determine the correct hypothesis, 
which is equivalent to determining the position of an edge in a column. 

Three different range estimates are obtained according to the three hypotheses. 
The estimate which matches the measured range best and is within the specified 
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error statistics should belong to the correct hypothesis. 
Decision Rule Mechanics 


If branch 1, Fig. 36, is found to give the best fit to the data, then it 
is chosen which implies that the edge occurs between stages i and i+1 in 
the column under consideration. Another decision tree is then initiated 
at stage i+3 to start reaching for the next edge, if any. If branch 1 
cannot be chosen, then it is rejected and the next observation is processed. 
At this next stage, previous hypothesis 2 (branch 2) becomes the new hypo- 
thesis 1 (branch 1) and a new hypothesis 2 (branch 2) is generated at stage 
i+3. This process is repeated till the end of the column is reached, f 

To work out the details of the scheme, the following two things are neces- 
sary. 

1) Estimates of ranges and slopes conditioned on the three 
hypotheses at each stage of observation. 

2) A decision rule which will determine the correct hypothesis 
after having obtained the estimates. 

Examining Fig. 37 it is obvious that the horizontal components of the edges 
(top and front edges of boulders) will be detected while processing columns. 
However, the vertical components of the edges (side edges of boulders) can 
be determined only if we process rows also, in the same way as columns. 
Therefore, a similar analysis is performed on the rows, Fig. 38. 

Model Definition 


The rangefinder scans an area in front of it by varying azimuth and the 
elevation angles to obtain an (m x t>) range matrix which has the same noise 
in it due to the errors in measurement, Eq. (1), Fjg. 35. The noise is 
assumed to be zero mean, Gaussian wi‘h a standard deviation In Eq. (1) 

Zjj is the noise contaminated measurement. 
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Model Along a Column 


( 2 ) 


Even though the adjacent rows and columns *re related due to their proximity. 











they may be considered independent of each other. The range (dj) and 
the changes in the ranges (g£> are defined to be the states of the system 
for each element of a single column. For a fairly gradual terrain g£ will 
be slowly varying, if not constant. With "i" the variable index along 
the length of any colunn, the system equation can be written as follows: 
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where 
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the quantity u^ is 


the change in g£ which will occur at the bottom edge of a boulder or the 
far edge of a crater. The Kronecker delta emphasir.es the fact that 

the change occurs at the edge at stage K. The quantity f . which is the 
ratio of gj + i to gj depends on the variation of ga along the terrain and 
f j > 1 if the column is processed upward, i.c. with increasing range on a 
flat terrain. It may be set equal to unity as an approximation. is a 
noise vector representing the roughness of the terrain and the inexactness 
of the model. 

On a flat terrain, the range, dj , will steadily increase, Fig. 36 and gj 
is almost constant. As the range readings cross the bottom edge, d^ will 
not Increase as much so that there will be a sudden decrease in g^. This 
is the abrupt change, i.e. jump in g£ to be detected. The sudden increase 
in d^ which occurs at the top edge of a boulder or the near edge of a crater 
can be modelled as two consecutive jumps of opposite signs in g^ and g^ + j 
Fig. 37. Eq. (2) along a column reduces to the following 

Z i+1 " ^*1+1 + ^l+i <4) 


where the error ‘IT'i+i has the statistics, S • 0 and VC ( v-^ V.) » 

q-'* £ lj. The matrix - ■ jl ojsincc only range measurements are obtained. 

Model Along a Row 

Only the horizontal components of the obstacle edge shown in Fig. 35 will 
be detected when columns are processed separately. The vertical components 
will be detected only if the rows are processed too. The vertical components 
in a row cause changes only in the ranges, hence it is unnecessary to estimate 
the range differences, and only a one-dimens tonal state is defln'1. 

Vi* W“it s jK J • *■" (5) 

is the change that occurs at the vertical component of an edge in any 
row. The covariance q is given some value after experimentation. 


Estimation of d« and in a Column 

If there are no obstacles In the path, the term Jg in Eqs. (3) and (5) is 


ia zero and Kalman Filter will give good estimates of the states. 

Occurrence of an edge will change d K or g K suddenly and KF will not 
then be optimal. In such case, an estimate of the change conditioned on 
the three hypotheses defined earlier will have to be calculated. 


Input Estimates 

For the two dimensional case, Eqs. (3) and (4) show that u t will appear 
in tj + 2 but not in z* or Therefore, a second residue is defined as 
follows: 


i+2 


, A 

■ z _ - H F F -V 
i+2 ^ i+1 i ^i 


( 6 ) 


and a suitiole performance index may be defined and minimized to show that 
r ^ + 2 is the best estimate of u^. 
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Similarly, u 1+ ^ 



K 


i+3 


( 8 ) 


New state estimates arc obtained by augmenting the KV estimates, ,y the 
input estimates obtained above 


X1 i + i ■ E<x m /H i> - V, * »“! 


^1+2 * E(X 1+2/B 2 ' * 1+2 + *Vl 

A A 

where and * r< KF estimates. 


(9) 


Top Edges of Boulders and Front Edges of Craters 

The parameters B and u^ in Eq. (3) are defined to model the changes in 
directly, which occur due to slope changes in the terrain. However, the 
top edges of boulders and front edges of craters cause large changes in the 
ranges "di" themselves. A Jisap ci d i4 » to d*.^ in Fig. 37 can be equivalently 
modelled as two consecutive changes of opposite polarity, l.e. to g 4 . 
•nd g f to g . g. These junps are always large as compared to the small 
Jmsps caused by the slope variations. Therefore, u^+j is also calculated if 
u* is found to be large. The program uses 2 meters as the threshold value, 
which, if exceeded by u^, indicates the occurence of a range jtmp. 


The Decision Algorithm 

Bayes decision rule. Ref. 11, is used to determi te which of the three hypotheses 
is most likely to be correct. The complete det .iis are described in Refs. 7 
and 8. 
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The horizontal component of the edge detected in a column is denoted by 
the symbol "x", the vertical component of the edge detected while moving 
along a row is denoted by "0" and their overlaps by "x" as in Pigs. 39-44. 


Results 


Figs. 39-44 give the outlines of boulders and craters at different dis- 
tances, with the scan concentrated around the obstacles. Most of the 
bottom edges of the boulders and some of the . ir edges of the craters are 
obtained. The vertical components of the edges denoted by zeros in the 
figures for the craters give an idea of the extent of the craters, as 
opposed to only the front edge obtained by th° FDR. 

Larger increments in the distances, a*- greater ranges, cause the bottom 
edge to be detected more clearly for the boulder at 30m than at 10 m. 


Conclusions 


By this method of added complexity, more information lias been extracted 
from the range matrix. The FDR computer program required 3 seconds to 
process the range matrix and obtain the incomplete edge wh*le the proposed 
RES method requires 6 seconds to obtain a more complete outline, on the 
IBM 360/67 computer. A choice or a compromise will have to be made depend- 
ing on (a) time between scans, (b) available computational capacity on the 
rover, (3) satisfaction of requirements, etc. 


Task D. l.b. Minicomputer Simulation of Obstacle Edge Detection Systems - 
J. Sher and D. Jacobs 
Faculty Advisor: Prof. C. N. Shcn 


The objective of this task was to develop a simulation of the laser range 
finder data of multiple obstacles on a minicomputer and to display visually 
the application of the methods of perceiving the hazard by the edge detection 
methods. This extends prior work which was directed to detection of a single 
hemispherica' boulder, kef. 12. 

Range Calculation 

Figure 45 shows *:he parameters ol the laser scanning range. Tie laser range 
finder is mounted on a mast, d meters above the ground. An on-board computer 
will control the direction of the beam, varying the elevation angle, P , 
and the azimuth angle, 0 , in order to form a searching area on the ground 
in front of the Rover. The total azimuth and elevation angle ranges are 
each divided into 25 equal increments. Starting at the minimus azimuth angle 
and at the minimus elevation angle, the range data is determined for that set 
of angular coordinates. The elevation angle it row increased by one increment, 
and the range data is determined again. When the range data has been detemiued 
for all of the elevation angles at the mlnimun azimuth angle, the azimuth 
angle is increased by one increment and the above procedure is repeated. The 
range data is stored in matrix fern, where each element of the matrix ropre- 
sentj the spherical coordinates of that point. 

The value of the v ange data for each pair of coordinates is equal to the 
distance along the laser beam to where the laser beam Impinges on the oust 'He. 
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Fig* 40. Outline of a Hemispherical Boulder of One Meter 
Radius at a Distance of Twenty Meters* 
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Fig. 41. Outline of a Hemispherical Boulder of One Meter 
Radius at a Distance of Thirty Meters 
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F.g. 42. Outline of a Hemispherical Crater of One Meter 
Radius at a Distance of Ten Meters 
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Fig. 43. Outline of a Hemispherical Crater of One Meter 
Radius at a Distance of IVenty Meters 
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Fig. 44. Outline of a Hemispherical Crater of One Meter 
Radius at a Distance of Thirty Meters 
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To find this distance, the height of the laser beam, starting at a minimum 
distance along the laser beam, is compared to the height of the obstacle 
at this point. If the difference between these heights is less than a 
given tolerance value, the length of the laser beam is calculated. If the 
difference is greater than the tolerance, the distance along the laser 
beam is increased by a given increment value. # The above procedure is re- 
peated until the difference is within the tolerance value, or until the 
sign of the difference changes. A sign change indicates that the obstacle 
is between the present coordinates and the one previous to the increment. 

A bisection search is now used to converge in on the obstacle - laser beam 
intersection point. 


Figures 46 and 47 show this searching scheme. Points a, b, c, etc. show 
the projection of the laser beam as it increments out toward the obstacle. 
Due to the long execution time of the minicomputer, the searching range was 
shortened from the minimum and maximum Y values of the laser to the minimum 
and maximum Y values around the obstacle, as shown in Figure 48. 


To simulate the occurence of noise, as would be found in the actual system, 
a range data matrix corrupted by noise is calculated. An algorithm des- 
cribed by Ref. 13 has been applied to generate a gaussian distribution of 
random numbers, which is then added to the range data matrix. 

Edge Enhancement 


In order to determine the coordinates of an obstacle from the range data 
that has been corrupted with noise, the edge of the obstacle is calculated 
using the normalized Laplacian algorithm, Ref. 14. This consists of com- 
paring an element in the corrupted range data matrix with the surrounding 
elements. The obstacle edge will be located at the coordinates where there 
is a large change in the value of an element in the rage data matrix. This 
sharp change of values is shown up using the following equation: 




1,J “i-l.J + “i+l,j +a l, j-1 + a l,j+l 


where a is an element of the range data. When the value of the magnitude 
1» J 

of the Laplacian is greater than a given threshold value, the obstacle edge 
is located at those coordinates. A matrix is now formed that indicates at 
what location the edge is at. 


The IDIIOM (Information Displays, Incorporated, Input-Out Machine) provides 
an interactive graphic display system. The parameters are displayed and the 
top view of the rover, obstacles, and searching area are displayed showing 
the relationships at these parameters. The program can be Interrupted by 
pointing the light pen at an area of the screen designated for the specific 
action to be taken. For example, a parameter can be changed by pointing 
the light pen at the area of the screen that has the value of the parameter 
displayed. This will cause a message to be printed on the typewriter 
terminal telling how to enter the new value for the parameter. 
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Results of the Simula t Ion 


Figure 49 illustrates the layout of the obstacles used in the simulation. 

The obstacles include a pyramid, a boulder, and a crater. 

The first simulation places the pyramid at (1.2, 17), according to the 
coordinate system in Figure 49. The boulder is at (-1, 20), and the crater 
at (0, 13). The height of the pyramid and the radius of the boulder are 
all 1 meter, the radius of the crater is 1.5 meters. The laser is placed 
at the origin with a height of d = 2 meters. The total azimuth angle is 
15 degrees, assuming that (©min ) and (©„, ax ) have the same magnitudes. The 
vertical scan covers the region between 10 meters and 50 meters from the 
; rover, corresponding to (B min ) = 2.3 degrees and (B,^) * 11.3 degrees. 

Simulation results from the IBM 360/67 are shown in Figures 50 and 51. 

Figure 50 includes the undistorted range data matrix and the one corrupted 
by noise. Figure 51 shows the normalized Laplacian matrix and the edge 
enhancement results. The threshold values used in this case are .03 for 
the inner edge and -.04 for the outer edge. In Figure 51, the inner edge 
is represented by "X", and the outer edge by 

In the second simulation, the pyramid is placed at (1,20) with a height of 
1 meter, the boulder at (0,11) with a radius of 1 meter, and the crater at 
(0,15) with a radius of 2 meters. All other parameters are the same as those 
used in the first simulation. Figure 52 shows the undistorted and the noise 
corrupted range data matrices. The normalized Laplacian matrix and the edge 
enhancement results are included in Figure 53. 

The simulation on the IDII0M has the advantage of a direct access to change 
the parameters. Three pictures can be displayed on the screen, as shown 
in Figures 54 through 56. Figure 54 lists the specifications for all the 
parameters used. All the parameters in this figure can be changed by using 
the light pen for interrupting and the typewriter terminal for inputting 
the new values. 

Figure 55 displays the layout of the obstacles and the terrain, correspond- 
ing to the specifications in Figure 54. Figure 55 shows the edge enhancement 
results. The second line includes the inner and the outer threshold values, 
which can be changed by using the light pen and the typewriter terminal. The 
bottom line is used to transfer pictures, e.g., by pointing the light pen 
at "top view", the top view in figure will replace the existing picture. 

By using the graphic display, the simulation becomes interactive with the 
operator. The parameters are displayed on the screen, and their values 
can be changed by interrupting the display with the light pen and typing in 
the new values on the typewriter terminal. Light pen interrupts can also be 
used to display the top view, display the range data, display the edge 
enhancement, display the parameters, or reset the parameters back to their 
initial values. 

This simulation provides investigators with a flexible tool which can be 
used to evaluate alternative schemes for detecting obstacles and the effect 
of data density and accuracy. 






Range Matrices R and R' of the First Simulation 
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Fig. Si. Matrix N and Edge Enhancement of the First Simulation 
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Fig. 52. Measurement Metrics R and R' of t ie Second Simulation 



Fig. 53. Matrix N and Edge Enhancement of the Second Slnulatlon 
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RECOGNITION OF 3-D OBSTACLES 

*** A SIMULATION 

> 

STUDENT : JONE-SU1AN SiiEU 
ADVISOR: PROF. C. N. SHEW 

THIS PROGRAM SIMULATES THE EDGE DETECTION SCHEME 
USED BY THE RANGE FINDER OF THE MARS' ROVER 


A PYRAMID, A BOULDER AND A CRATER ARE USED AS OBSTACLES 
OBSTACLE SPECIFICATIONS 


(I 1. LOCATION ( X,Y ) IN METERS 

I .' 

RO v SR IS AT ( 0 .0, 0 .0) 

[‘ PYRAMID IS AT ( 1 .0 , 16 .0) 

i ! B UDDER IS AT (-1 .0, 18 .0) 

CRATER IS AT ( 0.0, 12.0) 

I 2 . SIZE IN METERS 


* j. HEIGHT OF THE RANGE FINDER IS 2.0 M 

f ; HEIGHT OF THE PYRAMID IS 1.0 M 

! 1 * RADIUS OF THE BOULDER IS 1.0 U 

i RADIUS OF THE CRATER IS 1.5 M 

! i; 

.< 1-' . SCANNING RANGE AND HOIS”: 



1. HORIZONTAL SCAN COVERS 15.0 DEGREE 

MAX. THExh = 15.0 DEGREE MIN . THETA « -15-0 DEGREE 

2 . VERTICAL SCAN COVERS 10.0 M TO 50.0 M FROM ROVER 
MAX. BETA = 11.3 DEGREE MIN. BETA = 2.3 DEGREE 

I 

I 

3. GAUSSIAN NOISE: MEAN VALUE 0.0 CM 

STANDARD DEVIATION 5.0 CM 


RESET TOP VI EW RANGE L*.TA EDGE DETECTION 


Fig. 5A. Specification Display 
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RVR: ROVER 
P: PYRAMID 
Ik BOULDER 
C: CRATER 
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Task D.2. Obstacle Classification - K. Leung 

Faculty Adviser: Prof. C. N. Shen - — 

The objective of this task is to develop the minimum set of parameters and 
the conditions required to determine the passability or impassability of 
terrain. The ability to classify local terrain features by a minimum 
number of parameters is of critical importance to those groups concerned 
with the development of path selection algorithms, Task C, and of hazard 
detection hardware, Tasks A and B. 

It has been determined that a minimum of four parameters as related by 
three inequality conditions which must be satisfied simultaneously will 
be necessary. Further analysis may reveal that additional information may 
be required to achieve path selection systems of high reliability while at 
the same time minimizing built-in safety margins, i.e. as through unneces- 
sarily large threshold values. 


The analysis was approached by considering a number of fundamental questions. 
Is the object under consideration positive (boulder) or negative (crater)? 

Is the object of a size relative to the vehicle so that it deserves con- 
sideration, Figure 57? Are the heights and gradients of the object either 
too large or of the wrong shape, Figure 58? Is the object in a surrounding 
as to make it an obstacle, Figure 59? Does the combination of objects 
represent a hazard, even though the objects taken separately are not a 
hazard. Figure 60? 


From these and other related questions were deduced a set of Conditions of 
* Impassability. 

1. Instability - This is the situation where each wheel of the 
rover is subjected to terrain of different height such that 
the rover might be overturned. 


3 
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Shown in Figure 61 is a parameter CX , which is defined as 
the angle made by the rover's direction of forward movement 
with the horizontal plane. For a vehicle, there is a threshold 
value of C* , denoted as such that the vehicle is stable 

when c.< is less than this value. The value of * depends on 
the position of the vehicle's center of gravity and the dis- 
tance between the front and rear wheel axles, which can be 
predetermined experimentally. Hence, a prospective passable 
path must produce an less than along the path. 

Similarly, an angle <, can be defined as the angle made by 
the wheel axle with the horizontal plane and its maximum 

value for which the vehicle is stable, as shown in Figure 61. 

2. Steepness - The vehicle's limited traction will not permit it^ 
to climb a slope of more than r * 2 . The threshold angle 'X \ 
is due to stability consideration only. Since —2 is smaller 
than ^* 1 *. stability can be insured if “X?* is used as the 
threshold value of o< for passability. Hereafter, c< 2 *wiil 
be referred to as x'*. 

3. Abrupt Changes in Height - Because of the structure of the 
vehicle it cannot climb a vertical height of more than H*. 


► 








Fig. 60. This Combination of Boulder and Crater May Cause Instability 
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Hence a parameter H is defined as the vertical height 
of a step. Therefore if a prospective path is passable, 
the H along the path must be less than H*. 

4. Clearance - This is the problem that the underside of 
the rover may contact the terrain, either damaging or 
iraraobolizing the rover. 

Consider first a two-dimensional case. Assume that the 
dimensions of the wheels and their positions relative to 
the rover body are as shown in Figure 62. The distance 
between the wheel axles is denoted by L. Let S be the 
one-dimensioral terrain, given by 


z - f x (X) (1) 

Let the line segment L representing the wheel contact 
points be moved such that its end points coincide with 
the surface of the terrain. For points on S between P 
and Q, the distance between terrain surface and the 
segment L can be calculated. Depending on the design of 
the rover, it is necessary that adequate clearance be 
provided. If D is defined as the available clearance, then 
g < D for all points on the path trajectory. 

In the more general two-dimensional terrain problem, assur- 
ing that adequate ground clearance will be available is more 
difficult. It is necessary to fit a plane to the four wheel 
contact points and compare the terrain coordinates for all 
points within the space to this plane, The maximum value 
of the difference, g, so calculated is to be compared with 
'"he clearance, D, determined from the vehicle design. 


t 


t 





L 

i: 

u 

o 

e 


Lists of Situations, Parameters and Passabllity Requirements 

Below is a list of situations the rover will encounter as possible 
obstacles. For each situation ther-' is a set of parameters in terms of 
which the passability requirements 11 be written. 

S. : Vehicle is roving on flat inclined plane as shown in 

Figure 61. 

Parameters; Cn’ and v 

Passable if: . tan *.2 , tan > ' v 2 . 

< t*nK ,) + ^ 1 

2 o 

The sum of squares (tan v) + (tan v)* is the square of the gradient 
of the terrain the vehicle is on, which leads us to the above expres- 
sion. 


Sjt Vehicle is roving on boulders (craters) on a 
horizontal plane, as shown in Figure 63. 


x 


* arc sin 



► 


Parameters : 


w 




Fig. 63. Boulder (Crater) Combinations 




Fig. 65. Vehicle le descending a step 
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Passable if: (~;*) 2 + (^T^) 2 < 

t &n ^ tan , 

where is the pitch angle of the vehicle. 


1 


S^: Vehicle is descending a step, as shown in Figure 65. 

Parameters: t< - - arc sin q 

Passable if: (—S*'*) 2 + (— ^) 2 <. i 

tanct tanj 

where ^ is the roll angle of the vehicle. 

S^: Vehicle is ascending a step, as shown in Figure 64. 

Parameter : H 

Passable if: (jj*) 2 < 1 


S,. : Clearance as shown in Figure 62 

Parameter g 

Passable if : (&*) 2 ’Cl . g* * D 

8 

Now, let us define a situation S to encompass S^, 


parameters : 

«• * 

H and g and 

passable if: 

can* 2 
v tan* 9 

+ <“r *> 2 < 1 

tan 

and 

*»* < i 


a- i 

<§,)’< i 



S 9 ,..., and Sj, with 

( 1 ) 

( 2 ) 

(3) 


The reason for having three inequalities is that -<■ and £ are dependent 
on each other, H is Independent of c- , / , and g, and g is Independent of 
the other three parameters. The angle is the maximum value of only 
when y ■ 0, and y is the maximisa of y when *0. Hence the use of 
Inequality (1) eliminates the case when both >'and y are clc^e to but less 
than ex and respectively. 

Also, coupling Inequality (2) with inequality (3) using a logical AND * 
means^thot it is acceptable to have both H and g close to but less than H 
and g respectively. Therefore, situation S is not an obstacle only when 
Inequalities (1), (2) and (3) are all satisfied. 

As a next step, it remains to demonstrate that the inequalities which are 
proposed and the parameters required therein will provide the essential 
terrain information required for an efficient yet reliable hazard detection 
and path selection system. It may prove necessary to Increase the number 
of descriptive parameters and conditions if reduced safety margins are found 
to be necessary. 
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Task D.3. Terrain Modeling - Audbur Thompson 

Faculty Advisor: Prof. C. 'N. Shen 

In order for a Mars Roving Vehicle to navigate autoncmouslv , it must 
possess the capability of recognizing passable terrain, which in turn 
requires the ability to mathematically model that terrain. This task 
Involves the development of a method which will supply such a capability 
to the rover. The major task is to develop a scheme which cat' model 
small sections of the Martian terrain using laser range findf,r data. 

The scheme must also be able to limit discontinuities present at common 
boundaries of modeled sections. The purpose of such a procedure is to 
determine gross slopes and gradients of the nearby terrain. With this 
information, it will be possible to define lines ‘of constant gradient, 
which will identify forbidden zones where the gradient is too large for 
the rover to negotiate. 

Note that obstacle detection is not an objective of this task. The 
objective is to find general gradients. In fact, one aim is to find the 
gradient of terrain upon which small objects about the size of the rover 
are located. The detection of these objects if left to Tasks D. l.a and 
D.l.b. 

Concept of the Two-Dimensional Spline Function 

From a single scan a laser range fi.naer will provide a large number of 
data points to the rover. These points must somehow be fitted to a two- 
dimensional polynomial which models the terrain. One possible approach 
is to fit a high order polynomial expanded about one point to a large 
number of data points. But nigh order polynomials involve large numbers 
of cross product terms, and expansion about one point leads to less local 
accuracy. Thus the approach chosen for this task is the two-dimensional 
spline function. Sections of the terrain are fitted piecewise to low 
order two di lensional polynomials, each expanded about a different point. 
The terrain is then represented by a large number of such sectional poly- 
nomials. 

Initial Modeling Scheme 

In the initial analysis, two sections which share a common vertical 
boundary arc modeled. To formulate the spline functions, six data points 
arranged aj in Figure 66 are used. Each data point is given a cross-path 
coordinate, an in-pathcc .dinate, and a height, denoted by a, b, and h, 
respectively. The dashed vertical line represents the connon boundary. 

Consider the four leftmost points. * The height of any point within the 
area they enclose may be represented by the two-dimensional polynomial 

h - C o0 + C lfi (a-a ) + C 0 ,(b-b )+C,,(a-a )(b-b ) 

In finding the C's, which are unknown coefficients, the given heights of 
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Fig. 66. Data Point Arrangement for Initial Modeling Method 
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Fig, 67a. One Section 
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Addition of a Fourth Section 
Having Two Common Boundaries 
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the 

four 

data 

points may be expressed in matrix form as 
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Next consider 

the 

rightmost four points in Figure 66. Since it is 

desired to e)q>and 

about (ajotbio^* 

the matrix equation for these 

point 

heights is 
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In order to proceed later with the optimization step, both systems must 
be expressed in one matrix equation. Consequently, 


h# ■ Ax where h# 


h’ 

j 


o i 

X -! x ’ 

lh* 

1 i° 


x | X* 


The height of the terrain anywhere within these six data points is 
represented by the pair of two dimensional polynomials. Along the common 
boundary, the height, cross-path slope, or in-path slope of any point ma; 
be determined using either section model. It is thus possible that dis- 
continuities exist at such points for these terrain parameters. 


Method of Optimization 


Discontinuities in slope values will make it difficult to obtain reasonably 
accurate iso-gradient lines, and necessitates the development of some method 
which tfill limit these discontinuities. The approach chosen is the minimiza- 
tion of a performance function comprised of three terms. The first term 
indicates the square of the difference between the measured height of the 
data poirts, and the height at the same points of any optimized model. The 
second and third terms are, respectively, a measure o* the discontinuity 


occurring in cross-path and in-path slope along the boundary, 
scalar weighting and p , the performance function becomes 


Including 
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where R and L refer to '.right and left, 
and expressed In matrix form so that 


The last two terms can be expanded 
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If J is minimized with respect to X, the new estimate of che coefficients, 
now optimized, is 


X 


A A + 


+ JOG 


1 A T h 


Stochastic Problem 


If range measurement noise is considered, a minimum variance calculation 
should be perfotmed in conjunction with the preceding operation. This in* 
voives the inversion of an 8x8 matrix expressing the covariance of the 
heights of the data points used. Bur. since both sections are dealt w. r>i 
simultaneously, ~~'i two points are in conanon to the sections, two rows of 
the covariance matrix are identical. Thus the matrix cannot be inverted. 
This stochastic analysis fails because eight unknowns arc being determined 
simultaneously from six pieces of data. Due to this impasse, it was 
decided to modify the scheme. 


The Revised. Extendable Model 


The revised modeling scheme and optimization procedure are simply modifi- 
cations of those described previously. The aim of tho sodding schetr 
to mai. a in the t >tal number of unknowns found to he less than or equal to 
the number of additional data points used. For the four term polynomial 
used previously, a square array of nine points arranged as in Figure 67a is 
used. As noted in Figure 67d, the minimus number of data points added in 
expanding to an adjacent section is four. For a given section, ths four 
unknowns can be found iron the nine data points by using the least squares 
method. The appropriate matrix equation for the heights is 
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The revised optimization procedure works with only one see of coefficients 
at a time. To do this, assume that the elements of vector X have been 
found, and arc fixed. Consider a data -action to the ismedlate right, as 
shewn In figure 67b. In a manner analagoua to the one above, a matrix T 20 
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and r vector h20 ma y be formed. The resulting performance function, 
which ^als only with the elements of X 20 is similar to the one used 
initially, except a term has been added to penalize boundary height dis- 
continuities. 



Remember; ag that the elements of X 00 are fixed, this expression may be 
written in matrix form, minir. ized with respect to X£g, an optimized set 
of coefficients. 


If a section directly above the firrt is considered, as in Figure 67a, 
a set of new, but similar boundary conditions may be derived. 

Further to find X 22 for the additional section depicted in Figure 67d, 
it is only necessary to derive a similar performance function consisting 
of six boundary terms. Six terms are needed since two common boundaries 
exist. Note that for the resulting optimization, Xq 2 and X 2 Q are assumed 
to be fixed. 


If it is assumed then that once found, the parameters of a section ex- 
pressed in a vector X are fixed, a recursive scheme exists which will add 
any number of additional sections to the model. Section polynomials can 
be found individually, without altering those of previously modeled sections. 


Numerical Example 


To demonstrate the modeling procedure, a numerical example which modeled 
six terrain sections was run on the computer. The data points are spaced 
one meter apart. The situation represents a planar slope on which is 
located a one meter radius boulder. The data point heights are as follows, 
where the dashed lines indicate common boundaries. 
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The best way to display the results is by the Iso-'gradient maps of Figure 68. 
For Figure 68a, the slope discontinuity penalty was given a zero weighting. 
The line made up of "2's" represents the gradient line whose value is 25°. 

The area Inside the two lines thus has a gradient greater than 25°, indicat- 
ing the presence of the boulder. \ 
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Fig. 68c. Slope Penalty » 10 
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But it is desired to find the gradient in the vicinity of the boulder. 
For Figures 68b and 68c the slope discontinuity penalty weighting is 
raised to unity, and then ten. For these higher weightings the 25° line 
disappears. Thus the model is producing a terrain representation whose 
gradients are those in the vicinity of the bpulder, which in reality is 
about 14°. 


Conclusions 


The proposed terrain modeling procedure can find general slopes and 
gradients of terrain, and subsequently define lines of constant gradient. 
The method is recursive, and any number of sections can be easily modeled. 
Each section requires only about .15 seconds of computer time on the IBM 
360/67. With further development it appears that the two dimensional 
spline function modeling technique can be used to formulate iso-gradient 
maps of the general terrain features. Beyond this point, the following 
extensions should be considered; 

a) A full stochastic study should be undertaken to 
determine the exact effects of range measurement noise. 

b) The present scheme requires rectangular arrangement of 
data points. Since this will seldom be the case for 
real terrain, the boundary conditions must be generalized. 

c) From other examples, it was found that higher slope 
penalties are needed in some situations than in others to 
attain adequate results. Thus it might be appropriate to 
study the possibility of variable penalty weightings. 


D.4.a. Measurement Scanning Concepts - Mark Friedman 

Faculty Adviser: Prof. C. N. Shen 

The objective of this subtask is to develop scanning concepts for 
laser rangefinder hazard detection system. This involves determining the 
minimum field of view and data density. This is required to ensure that 
the measured range data will provide accurate and reliable results when 
used by obstacle detection and the terrain modeling techniques described 
under Tasks D.l.a and D.3. 

The basic scanning parameters are shown in Figure 69. The field of view 
is defined by P , the rangefinder elevation angle measured down from 
the horizontal, and 0 , the rangefinder azimuth angle. The rangefinder 
will scan the field at constant increments of ,5 and 0 denoted by 
and <1© respectively. The inpath point spacing is &b and the cross- 
path spacing is .la. The quantity -Ab ' is the vertical data point 
spacing as shown in Figure 69. 

The safety requirements and capabilities of the rover provided the criteria 
for defining the field of view parameters given in Table 1. The require- 
ments of the terrain modeling and obstacle detection procedures made it 
necessary to divide the scan into a near scan (3- 16m) and a far scan 
(16-30m). Each was treated as a separate field of view and using the same 
criteria as above, the parameters of Table 2 were derived. 
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Table 2 

Divided Terrain Scan 



lieu** Scan 

Far Scan 
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57.5° 

10.6° 

wax 

6 mln 

-lk.C° 

-19.6° 

^tet 

72.1° 

30.2° 

*U* 

53.0° 

53.0° 

A0 

X.U° 

.76° 

AO 

3.6° 

1.91° 

N 0 

15 

28 

N 8 

52 

kO 

N 

780 

1120 


N * 1900 data points 
toi 
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The data point spacing parameters were determined from the requirements 
of the modeling procedures. However the obstacle detection routines 
call for a significantly smaller data point spacing than the terrain 
modeling procedure. This problem was approached in two ways. 

The first procedure referred to as the Terrain and Zoom Scan Procedure 
is to scan the field with a data point spacing appropriate for the terrain 
modeling technique. This spacing would not give a well defined outline 
of an obstacle but would at least give an indication of its presence. 

The second scan would be a "zoom" scan with reduced spacing aimed at the 
obstacle location indicated in the previous terrain scan. The total 
number of points in the terrain scan would be about 1900 (Table 2), Using 
an expected scanning rate of 1 msec /data point, this would give a scan 
time of 1.9 sec. The time of the zoom scan would probably be from 2 to 
5 sec depending upon the size of the obstacle. The disadvantages of this 
method are considerable. The hardware modifications for the zoom scan 
are likely to be quite complicated and the presence of more than one 
obstacle in the field would require multiple zoom scans between each terrain 
scan. 

These disadvantages may be easily avoided by again using near and far scans 
with data point spacing set for obstacle detection. Software techniques 
may then be used to systematically skip over scxne points to achieve the 
spacing required for terrain modeling. 

In an effort to reduce the total scan time, three sets of far scan para- 
meters have been derived for this obstacle scan procedure and are .given 
in Table 3. Each saving in scan time was made by reducing the field of 
view of the far scan. 

The scan time and total number of data points taken can be reduced by a 
third procedure namely the Corridor Scan dividing the field into M 
corridors as shown in Figure 73. This procedure is to be applied to the 
full width for ranges from 3 to 30 meters with M * 5. The scan will again 
be divided into two subscan types. First, an initial scan with low data 
point spacing will be taken of the whole field to merely indicate the 
presence of obstacles. Then a corridor scan with higher data point density 
will be taken for only those corridors which show premise of passabillty. 

The Corridor Scan method is expected to result in reduced necessary on- 
board computer core and reduced calculation time of the range data pro- 
cessing techniques. Tables 4a and 4b. 

By examining the avoidance time and the number of scans per field for the 
Corridor Scan Procedure, it is seen that the maximum scanning range can 
be reduced from 30 meters with little added danger to the vehicle if the 
vehicle's average speed is less than about 0.2 meters per second. 

Terrain and Zoom Scan Method 


The rover is expected to be able to traverse a positive or negative slope 
up to 25°. This is used to calculate ( 9 . The quantity 5 max * 8 
for the situation of the rover approaching a negative slope change of 2S° 
and scanning at 3 m for the near scan and IS meters for the far scan as 
shown in Figure 70. Figure 71 shows the situation use£ for the calculation 
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Fig. 70. Situation for Beta Max. 



Fig. 71. Situation for Beta Mi tv 
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of ^ nln . The rover is approaching a positive slope change of 25° 
and scanning at 15 meters for the ner scan and 30 meters for the far 
scan. 

The quantity 9 is calculated so that the rangefinder can scan the width 
of the vehicle, expected to be 8bout 2.1 meters plus an allowance on 
each side or 3 meters total at a distance of 3 meters as shown in 
Figure 72. Table 2 gives the value, derived above. The values for the 
field of view parameters in Table 1 and Table 3 are similarly derived. 

In Table 1, A 9 was set so that the rangefinder would sense an object of 
one meter width with at least one data point at the farthest range. The 
quantity was s 't so that an object at least 0.4 meter high could be 
detected with one a. -a point at the farthest range. The height 0.4 meter 
was chosen because the vehicle is expected to be able to roll over an 
object of lesser height. 

With the scanning increments calculated above the data point spacing 
becomes very small at close ranges. This would give rise to large error 
in the terrain modeling procedure. Ref. 15. Even when the field of view 
was divided as in Tabic 2 the error was unacceptable for ranges closer 
than about 8 meters. 

The number of data points per scan is denoted in Table 2 as N. The total 
number of scanning points per field divided by the expected scanning rate 
of 1 msec /data point gives the total scan time of 1.9 seconds. 

If a terrain scan with parameters similar to Table 2 were used, a terrain 
modeling procedure. Ref. 16 would of course be used but an obstacle detec- 
tion routine would also be applied to the range data. Although the data 
point spacing would be too large to give a well defined outline of an 
obstacle, the routine would give an indication of either the top edge of 
a boulder or the near edge of a crater. Ref. 17. A soon scan would then 
be utilized with reduced data point spacing aimed at the obstacle loca- 
tion given by the terrain scan. The extent of the "zoom" would depend on 
the size and range of the obstacle. The number and range of the obstacle's 
warning points could be used to determine the proper amount of data point 
reduction. 

O bstacle Scan Method 

Table 3 gives the scanning parameters of the Obstacle Scanning Procedure. 
The scanning increments ££a nd 9 are calculated by setting As and lb* 
(see Figure 69) both equal to .2 meters at the farthest range for the near 
ard for the far scan. This gives a minimum data point density on the face 
of an obstacle of 5 points per meter which is aestased to be sufficient for 
a well defined outline of the obstacle. 

Three sets of far scan parameters have been derived and given in Table 3 
in an effort to minimise the total scan time by reducing the field of view 
of the far scan. Far scan A extends from 15 to 30 swters with a full 
angular width of 53°. Far scan B also extends from 15 to 30 meters but 
the angular width has been reduced to about 37° giving a minimum linear 
width of 10 meters. - The angular width of far scan C is reduced as in far 
scan B but extends from 15 meters to only 25 meters. From the values given 
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for the number of scanning points It Is seen that the total scan time 
for this method will be from about 11 to 16 seconds. 


Corridor Scan Concept 

In the corridor scan concept the Initial scan will have a low 'entity 
spacing. Its purpose is merely to Indicate which corridors contain 
obstacles but not to show their outlines. It will do this by comparing 
the range of each data point to range values of clear level terrain vhich 
will be stored in the on-board computer. Each point which differs signi- 
ficantly from its "clear" value will be assigned a danger value dj* 
where i Indicates the i c ^ corridor and j runs from one to N, the 
nunber of points in the corridor. Then each corridor may be assigned a 
danger value D^: 


D 


i 




i - 1, M 


The d values should be dependent on J. , the level of the obstacle 
point height and R the obstacle point range, with a closer and taller 
obstacle point being assigned a larger danger value. 

Aa an example of how X nay be determined, consider Figure 74. Range 
values for three levels are stored and compared to the data. If the 
obstacle point falls between -0.2m and 40.2m it is given J »0 and does 
not contribute to d. If it falls between 0.2m and 0.4m it is assigned 
4*1 and if it is above 0.4m or below -0.2m it is given 2. A higher 
danger value is given to negative obstacles since their actual depth is 
usually difficult to sense. 

The next step In the procedure is to make a scan decision based on the 
corridors' danger values. One possible method is as follows: 

(1) If there is one or more corridors with D u D g 
travel in the direction of the safe corridor nearest 
the center of the field with no detail corridor scan. 

(2) If Dj V D u ; i«l,M, the vehicle may have to stop and 
do a complete detail scan. 

(3) Otherwise choose smallest D, and do a detail corridor 
scan. 

(4) If detail corridor scan shows an obstacle, try next 
larger D. 

In step 1, D a is a low danger value such that a corridor with a D less 
than D s would be almost completely clear and level. A corridor near the 
center of the field is preferred since more is known about the terrain 
surrounding this cori idor than an outside corridor such as corridor 1 or 
5 (Figure 73). An exception to this rule would be If three corridors 
together each had D leas than D s such as corridor 1, 2 and 3. In this 
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case corridor 2 would be the safest. 

In step 2, Du is a greater danger value such that tncre Is a high proba- 
bility that there exists an obstacle which could not be avoided. Two or 
possibly thiee corridor scans may be necessary (step 4); however, if the 
terrain is so rough that three corridor scans show unavoidable obstacles, 
step 2 would probably have been satisfied and this field would have bs 
completely avoided before any detail scans were made. 

As shown in Figure 75, each corridor is narrower than the vehicle's width 
of 2.1 mcte.*s at ranges closer than 11.4 meters, so that there will be 
some ground ov^r which the vehicle will travel which cannot be completely 
modeled with the corridor scan alone. This gives a criterion for deter- 
mining the proper density of the initial scan. The density should be 
set at 5 points per meter at the range of 11.4 meters so both terrain 
modeling and obstacle detection techniques may be performed on this near 
sector. This gives the initial scan parameters shown in Tabic 4a. 

The parameters of the corridor scan as shown in Table 4b were derived by 
setting the density at 5 points per meter at the 30 meter range. Each 
corridor scan need only extend from 30 meters down to 11.4 meters since 
for closer ranges the initial scan density is proper for the ter: jin 
modeling techniques. As before, it is assumed that software pre-processing 
will be used to systematically skip data points to achieve the proper lover 
spacing required by the terrain modeling procedures. 


Scan Time Performance 


The maximum scanning range should be chosen so that, once the vehicle gets 
within this range of an obstacle, it should have enough time to detect it 
reliably and avoid it safely. This analysis so far has assumed a maximum 
scan range R ^y of 30 meters. This is somewhat arbitrary since the 
vehicle's average speed has not been specified. 

To ensure accurate detection of an obstacle, it is desired to have the 
number of scans denoted as >, which can be made as the vehicle travels 
the length of its field of view as large as possible. This number 
can be simply calculated from the following relation: 

(maxim up scan range) 

£*• » (scan time) (vehicle speed) 

This calculation has been performed for the Corridor Scan Procedure. 

Table 5 shows the results of this calculation for a vehicle speed of 1 
meter per second. The maximum range is R^* . T is the total scar: time, 
and Z is the distance traveled by the vehicle during one total scan. The 
total scan time is equal to the constant initial rcan time of 4.1 seconds 
plus the time of the corridor scan which is dependent on the maximum scan- 
ning range. Figure 76 shows a gt*ph of 5k versus R^y for various 
vehicle speeds. 

For safe obstacle avoidance, the time from the end of the total scan 
which covers the obstacle to the time at which the vehicle reaches the 
obstacle should be as great as possible. Assuming the scan begins, just 
when the obstacle is at K ^y , this time, called T t , is just equal ro: 
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TABLE 5 


Scan Time Performance 


B 

max 

(meters) 

T 

(sec) 

Z 

(meters) 

Q 

t a 

(sec) 

15 

1*.7 

1*.7 

3.19 

10.3 

20 

5.1 

5.1 

3.93 

14.9 

25 

5.7 

5.7 

4.38 

19-3 

30 

6.6 

6.6 

4.51* 

23.4 

35 

7.8 

7.8 

. 4.48 

27.2 
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9.0 

9.0 

4.44 

31.0 

For Corridor Scan 

Procedure. 

Vehicle speed equals 

1 meter per sec 
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Table 5 also shows the avoidance times for each for a vehicle 

traveling at 1 meter per second. Figure 77 is a grapnof T^ versus 
®max for various vehicle speeds for the corridor scan procedure. 

In Figure 76, it is seen that Q , the number of scans per field is 
indeed a maximum for a maximum scan range of 30 meters. However, for 
vehicle speeds of 0.2 to 1 meter per second, this performance is not 
greatly reduced by reducing the maximum scanning range to 20 meters. 

But at these speeds, Figure 77 shows that the avoidance time would be 
cut down by about 35%. For slower speeds on the order of 0.1 meters per 
second, the avoidance time is not significant in this consideration. The 
number of scans per field and the avoidance time is so great for a vehicle 
speed as low as 0.1 meters per second that the maximum scanning range 
could probably be cut all the way down to 15 meters or below with little 
danger to the vehicle. 


Conclusion 


The total scan time for the terrain and zoom method will depend on the 
number and size of the obstacles in the field of view of the rangefinder. 

For three obstacles and assuming the scan time of the zoom scan to be 
the same as the terrain scan, the total scan time would be about 8 seconds 
which is substantially less than the scan time of the obstacle scan approach. 
However, the hardware and software problems involved in the variable spac- 
ing reduction of the zoom scan and the aiming of the zoom scan are expected 
to be quite great. It is felt that the obstacle scan method, where the 
field is scanned with small data point spacing for obstacles and data 
points are systematically skipped over prior to terrain modeling, will be 
much simpler and more desirable. 

The corridor scan method is a combination of the terrain and zoom scan 
method and the obstacle scan method and avoids some of the disadvantages 
of each. For the obstacle scan about 15,500 points are taken in a scan 
time of 15.5 seconds, but an initial scan and one corridor scan takes only 
6.6 seconds for the necessary 6600 points. The scan time is reduced to 4.1 
seconds if the terrain is completely clear and a corridor scan is not needed. 
Besides the reduction in scan time, the Corridor Scan procedure reduces the 
necessary core required in the on-board computer and the complicated data 
processing techniques are used less often and with less data point. The 
hardware changes are expected to be small since the parameters of the 
corridor scan are fixed and a continuously variable data point spacing is 
not necessary. 

Once a vehicle speed is specified, the maximum scan range may be chosen 
to best reflect safety considerations using the curves of Figures 76 and 
77. Perhaps the maximum scan range can be reduced from 30 meters so that 
less data points will be acquired. If the vehicle's average speed is only 
0.1 meter per second, the maximum scan range can probably be safely cut in 
half. Then less computer core will be required and the range data process- 
ing techniques will be applied to less points, further simplifying the 
terrain sensing system. ~ \ 



Task D.4.b. Laser Rangefinder Error - M. Friedman- 
Faculty Advisor: Prof. C. N. Shen 

It is the objective of this subtask to investigate the range error due 
to the finite width and spread of the rangefinder laser beam. This is 
necessary in order to ensure applicability and reliability of the terrain 
modeling and obstacle detection procedures which use the range data. 

The beam, originating from a height of 3 meters, will strike the surface 
of the ground at a fairly acute angle at the farther scanning ranges, 
causing the projection to be spread along the surface of the terrain. 

Parts of the beam will return to the rangefinder receiver at slightly 
different times since they were reflected at slightly different ranges. 
Hence, even if a very sharp pulse was transmitted, the received pulse 
will be spread out significantly in time giving an ambiguous range 
measurement. 

It was found that the spread at the projection increases as the target 
range increases. The spread for a particular range may be minimized by 
proper choice of the transmitter aperture diameter. For the greatest 
range of interest, 30 meters, the optimum aperture diameter was found to 
be .78 centimeter. Table 6 gives the spread for various ranges between 
3 and 30 meters. The aperture diameter used was 0.78 centimeter which is 
optimum only at 30 meters. 

The measured range value will always be less than the true value by an 
amount dependent on the range. However, this analysis was performed with 
the assumption at level ground, i.e., the target is assumed to be oriented 
horizontally with respect to the rover. Therefore, the error cannot simply 
be added to the measured value to get the true value since the rover will 
encounter targets at various orientations at all ranges. 

If a constant discriminator set at a level of 50% is used as the return 
pulse time detector, the range error at 30 meters due to the beam spread 
will be about 6.27 cm. This accuracy is expected to be acceptable for 
reliable use of the image enhancing procedures. 


The Spread Equation 

Figure 78 shows the laser rangefinder geometry with the rover on level 
ground. The beam, transmitted from the top of the rover mast, may be 
pictured as originating from a point source located behind the rover. The 
beam is assumed to be a right circular cone of solid angle 2 ©^ where ©^ 
is the angle of dispersion of the beam. The beam diameter at the trans- 
mission point at the tope of the mast is equal to the diameter, D, of the 
final aperture of the transmitter optical system whether it is a lens, a 
stop, or the laser element itself. The intersection of the beam cone with 
the plane of the ground defines the projection as an ellipse whose length 
in the major axis is taken to be, S, the spread of the beam projection. 


From this geometry it is seen that S is a function of the mast height H, 
the range of the target Rq, the aperture diameter D and the angle of dis- 
persion ©d. The spread equation is: 
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and making small angle approximations for 0^: 


4 9. R 
d o 


*w 


Intensity Pattern 


If the angle of dispersion can be given in terms of basis rangefinder 
system parameters such as the aperture diameter and beam wavelength, the 
spread and therefore the range error may be minimized by proper choice of 
these parameters. This necessary relation is dictated by the beam intensity 
pattern. 


Theintenslty at a general point in the pattern is a function of 0, the 
angle from the beam axis to this general point, Ref. 18 and 19, as shown 
in Figure 79. Specifically: 


K«> 


. OTp sine ) 
1 N 

TT D sine 


where is the first order Bessel function 
D is the aperture diameter 
and X is the beam wavelength 
Figure 80 shews a graph of this relation. 


The first minimum occurs when the argument of the Bessel function is equal 
to 1.22 TT . If the beam is approximated by a cone extending only to this 
first minimum, the angle of dispersion may be determined by setting the 
argument of the Bessel function equal to 1.22T[ . The angle will be on the 
order of milliradians; thus, slnOj may be approximated by 0^. Solving for 


1.22 X 


Diameter Optimization 


Substituting equation (4) into equation (1) gives: 
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The mast height has been set at 3 meters, the range is to vary from 3 to 
30 meters, and the wavelength will presumably be set according to the 
actual laser used and will probably be in the infrared region. This leaves 
only D, the aperture diameter, to be adjusted so as to minimize the spread. 
The optimum diameter can be found by setting the partial derivatives of S 
with respect to D equal to zero. The second derivative is positive in- 
dicating that the spread will be a minimum for this optimum D. 


Pulse Shape 

As stated previously, the projection will be an ellipse. Figure 81 shows 
the ellipse and defines the coordinate system to be used in the pulse shape 
determination. The equation of such an ellipse is: 

+ -j£- - 1 (6) 


where h, a, and b are shown in the flture. 

These elliptical parameters can be related directly to the system parameters. 

If this ellipse is divided into rectangular sections as shown in Figure 81, 
all the energy returned from a particular rectangle will arrive at the range- 
finder receiver at approximately the same time. Therefore, each rectangle of 
width Ax is related to a unique time at the receiver. The Intensity of 
the light at the receiver at each t is proportional to the power density of 
the light leaving each rectangle. The constant of proportionality is assumed 
to Include the path length even though each rectangle is at a slightly dif- 
ferent distance. At distances on the order of tens of meters these small 
path-length differences have a negligible effect on the light intensity. The 
major effect of these slight distance variations is the varying time of flight 
of the energy reflected by each rectangle. 


The power density leaving each rectangle is equal to the product of the 
lntenslt) of the light Incident to the rectangle and the area of the rectangle. 
However the Intensity is not constant on any particular rectangle. Thus, an 
Integral of the following form will have to be calculated for each rectangle: 
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where I R - Intensity at re 


P - power density leaving target 
t 

K - constant of proportionality 

x - width of rectangle 


(7) 


and the length of each rectangle, J , is taken to be equal to the width of 
the ellipse at the value of x corresponding to the center of that rectangle. 
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Prom equation 6: 

r <«d>l 2 ] 4 

9 » 2y - 2b [ 1- a 2 J (8) 

If Ax is made the same for every rectangle the integration over x can be 
included in the proportionality constant. 

I t (x»y) is indirectly given by equation (3). 1(0) can easily be trans- 

formed from the polar system based on the circular beam to the rectangular 
system based on the elliptical spot. 

The Intensity of each rectangle can new be related to a distinct time of 
flight At. If each Alg is evaluated and plotted against its respective 
A t, the resulting curve will approximate the shape of the return pulse. 

Figure 82 shows the return pulse shape for a range 30 meters and optimum 
aperture diameter of .78 cm. Equation 7 was used with Ax set at 1.5 ssa 
giving about 100 rectangles. 

Subtracting the time at which the leading edge of the pulse reaches 507. 
of the maximum from the peak time gives a time error of .209 nanosecond 
or an equivalent range error of 6.27 cm. 

The optimal laser rangefinder parameters have been determined for the 30 
meter range. As the studies under Tasks C and D proceed, alternative 
maximum ranges of interest may be identified and a re-optimization of para- 
meters will be in order. 


Task E. Chemical Analysis of Specimens 

One important phase of the initial missions to Mars is the search for 
organic matter and living organisms on the martian surface. The present concept 
for attaining this objective consists of subjecting samples of the atmosphere and 
surface material to certain chemical and biochemical reactions and thereafter 
analyzing the products produced, probably in a combination gas chromatograph/mass 
spectrometer (GC/MS). The gas chromatograph is proposed for separating complex 
mixtures evolved from the experiments into small greups of similar chemical species. 
Chemical analysis of these groups would be accomplished in the mass spectrometer. 

It is the objective of this task to provide engineering techniques and design 
criteria for creating and improving such a system. 

Most of the previous effort has Involved the systems analysis of the 
gas chromatograph using simulation. Ref. 20, 21, and 22. This technique uses 
mathematical models, which incorporate fundamental parameters evaluated from re- 
ported experiments, to explore various concepts and to direct further experimental 
research. Application of prior work to characterize chromatograph performance and 
improve the mathematical model was considered during the past year. Task problems 
mere attacked by a two-member team, each of whom pursued a specific asslgnaant: 

1* Chromatograph system characteristics 

2. Chromatograph simulation development 



Task E.1 


Chromatograph System Characteristics - C, J. Rose 
Faculty Advisor: Prof. P, K. Lashraet 
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It was the objective of this task to develop techniques for designing 
chromatograph systems using simulation models and concepts briefly 
considered previously. Ref. 23. As the purpose of the chromatograph is 
to separate chemical species, this initial effort involved a quantitative 
definition of the degree of separation attained in a chromatographic 
column of particular design, operating under specific conditions. A large 
value for this degree of separation implies a chromatogram having distinct 
peaks, whereas a low value suggests poor separation or overlapping peaks. 

The degree of separation is a function of the physical par* ers of the 
chromatograph and the chemical characteristics of the materials being 
separated. On this basis, this task considered the following: developing 

a generalised simulation procedure; determining the effect of design para- 
meters upon the degree of separation through simulation; and obtaining 
realistic characteristics of the sample injection pulse. Progress, de- 
tailed in a recent report. Ref. 24, is summarised below. 

To evaluate the concept of separation, a computer procedure was developed 
to simulate the chromatograph. Simulation was accomplished by numerically 
convolving a previously developed model, Ref. 20 and 25, with typical 
sample injection pulses represented by data. The generalized simulation 
procedure treated a maximum of five components so practical mixtures could 
be considered, especially in the region of overlapping peaks. The simula- 
tion used actual sample sizes, carrier gas compositions, and real and 
dimensionless time so the results could be directly related to actually 
observed system operation. 

In a chromatogram, the output signal is the sum of composition contributions 
of all components. Based on this concept, chromatograph performance, was 
defined as 

S t - A^/A O) 

where 

S. ■ separation of component i for the peak 
where the component predominates. 

A. ■ contribution of component 1 to the area 

under the peak where the component predominates 

A » total area under the peak of interest. 

In these studies, the areas were computed from .simulations _over a finite 
peak width of four standard deviations (4 ' f.u> (2) where, w (2) was the 
second time moment about the mean) centered about the peak mean. Prior 
studies have shown that at least 94Z of the chromatogram area lay within 
this region. Sj was the effective composition of component 1 in the 
peak where component i should appear. It ranged from unity for complete 
separation to zero for essentially no separation; l.e., the sample was small 
compared to amounts of other interfering chemicals. Thus each component was 
assigned a velue of S’ indicative of its separation from the other components. 





For excellent separation, a high value of S, 
would be expected. 


srhep* 0.85 or higher. 


This index of performance required analysis of system simulations. To 
develop preliminary estimates of S prior to such simulations, the appro- 
ximation of the output pulses by Gaussian curves was investigated. Since 
Gausslan-type curves are defined by the first and second txmj moments 
which are known from the characteristics of the input pulse and the chroma- 
tograph differential equations, Ref. 25 and 26, the separation S can be 
rapidly estimated from design parameters without a simulation. 

For a binary system, assigning superposition of the individual chromatograms, 
the separation of component I became 1 


Aj + A n 


where the chromatogram areas A* and Ajj were obtained by simulation. The 
time interval for determining the areas was 


®I ± 2 V .a< 2 > 


where 0j was the mean or first time-moment about the origin for the chroma* 
togram of component 1 and pj(2) was the variance or second time-moment 
about the mean Op Rearrangement of the separation equation gave 

V A u - V<‘- S i> 


or as an approximation 


V A I1 


(Gaussian Area ^ /Gauss ian Area^) (Nj/Njj) 


where N. and N_, were the total amounts of components I and II res- 
pective^. 

For a given binary system under consideration, the first and second moments 
for the chromatograms for both components were determined as a function of 
the physical parameters from the system equations. Ref. 25. As shown in 
Figure 83, the Gaussian area for component 1 became 

A + ? 'i 


f (8) d 0 
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where 

0 ■ dimensionless time - vt/L 

L ■ column length 

t ■ time 

v * Interstitial gas velocity 


and y_(0) — was a normal distribution function having mean ©_ and 
variance Uj(2). For j * 2, the area was 0.9545. As shown in Fig. 83, 
the Gaussian area for component II included within the above integration 
interval was 


Gaussian Area^ ■ 


A + \ v^i< 2 > 


y n (e) d o 


J °i * \ ^ 7i<*> 


„ + \i w if n< 2) 

y n (6) d 0 
' * * n (2) 


where yn(6) was a normal distribution function having mean Ojj and 
variance Ju jj( 2). The time corresponding to ,Aj of component I was 
associated with a ?ijj(l) of component II. Similarly -^i was associated 
with a ^ii(2). For the situation of Fig. 83, the governing relations were 


- [»I - « n + V 7 j( 2) J / y^ u (2) 

>n (2> - [»j - 5 It - \ i ^ t ( 2 ) / 


Figure 84 has compared the approximate method of computing separation, 

Eq. 2, with the exact method, Eqn. 1. Simulations reported earlier, Ref. 23, 
were the basis of this comparison. Use of Gaussian-type outputs yielded a 
first-order approximation to the separation, although on the average the 
predictions were optimistic. However, it is a valuable preliminary design 
tool. 

To relate more clearly this research to the Mars mission, information was 
compiled on physical characteristics of proposed systems and possible 


Separation from simulation 



. Separation assuming Gaussian - like chromatograms 




Fig. 84. Comparison of separation functions for binary 
system derived from simulation and Gaussian - 
like chromatograms. 

0 - separation for first - appearing component 
A - separation for second - appearing component 
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chemicals postulated for separation. Research performed at Jet Propulsion 
Laboratory provided a guide on chemicals which might be expected. The 
chromatogram for a 24-component test mixture using a Chromosorb-W column 
coated with 5% Dexsil, Ref. 27, has been reproduced in Fig. 85. Estimates 
of the thermodynamic parameter hiRq for the components have been summarized 
in Table 7. These provide a practical range of values of the parameter for 
future system studies. 

Descriptions of physical systems have appeared in several reports. Ref. 28, 
29, 30, 31. For these representative systems, the Peclet numbers, or dis- 
persion parameters, were estimated using a parameter estimating program 
developed earlier, Ref. 21. As summarized in Table 7, Peclet numbers en- 
countered in practical chromatographs range from several hundred for mini- 
columns to 30,000 for very long columns. 

Previous studies, Ref. 23, have shown that the sharpness of the sample input 
pulse often determines whether or not separation between two similar sub- 
stances will occur. To determine realistic characteristics of input pulses, 
experiments were conducted on the test facility, Ref. 32, to obtain the 
effect of sample size, carrier gas flow rate, and other parameters upon the 
input pulses. Samples of gaseous ethylene were injected into a composition 
detector using a minivalve, shown schematically in Fig. 86. Typical ethylene 
pulses obtained for three helium carrier gas flowrates are shown in Fig. 87. 
Moment analysis of pulses from sixteen experiments have been summarized in 
Figures 88 and 89. The peak mean and standard deviation have been reported 
as dimensionless ordinates through use of T , the time ethylene remained in 
the sample loop during the injection process. The abscissa is the ratio of 
sample size to carrier gas flow rate. The ratio of carrier gas pressure to 
sample loop pressure has been shown parametrically. These curves summarize 
approximate effects of des ,n parameters on the input pulse characteristics 
and provide a means of assessing the effect of a design upon the separation 
occurring within the chromatograph. 


Task E.2. Chromatograph Simulation Development - M.P. Feinberg 

Faculty Advisor: Prof. P. K. Lashmet 

It was the objective of this task to develop and evaluate a practical, 
nonlinear adsorption model to account for composition-dependent effects 
observed in previous research. Ref. 33 and 34. A system of three coupled 
partial differential equations was derived earlier to describe, these effects. 
Ref. 35. Excessive computation time due to the complexity of the system of 
equations, Ref. 35 and 36, indicated a simplification of the model was re- 
quired. Voytus, Ref. 25, derived a simplified model which considered axial 
interparticle diffusion, convection, and equilibrium adsorption/desorption 
using a linear isotherm. Assumptions for this model, designated as the equi- 
librium adsorption model, included negligible mass transport effects between 
the carrier gas and the adsorbent and no intraparticle diffusion effects. The 
resulting partial differential equation was solved analytically using a con- 
volution integral. Ref. 20. Although the model simulated many systems well, 
in a few Instances model predictions deviated consistently from actual data. 
The variation of the thermodynamic parameter ibRq in such systems with 
composition supported the idea that a nonlinear, composition-dependent iso- 
therm should be considered. Preliminary work on solving the resulting 
nonlinear equations begun earlier, Ref. 37, was continued. Activities, 
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Table 7 ♦ Range of Parameters for Typical Chromatographic Systems 
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RANGE OF THE THERMODYNAMIC PARAMETER mRo 
(JPL Multicomponent Test Mixture) 


Breadboard 


Compound 

Time (min) 

Dimensionless Time 

mRo 

Air 

1-1*5 

2.58 

0.632 

Furan 

2.92 

5.19 

0.238 

Methyl Furan 

U.93 

8.77 

0.1286 

Benzene 

7.1*1 

13.19 

0.0820 

2,5-Dimethyl Furan 

9-53 

16.96 

0.0626 

Octane 

15.67 

27.90 

0.0372 

Ethylbenzene 

39-25 

31*. 27 

0.0301 

Nonane 

19-98 

35-57 

0.0289 

Decane & Mesitylene 

23 . 1*0 

1*1.66 

0.02h6 

o-Cresol 

25.00 

1*1*. 51 

0.0230 

Undecane 

25.68 

1*5.72 

0.0221* 

Dimethylanal ine 

26.70 

l*7.5l* 

0.0215 

Benzyl Cyanide 

29.17 

51-93 

0.0196 

Tridecane 

30 . 7 ' 

5l* 66 

0.0186 

Indole 

32.30 

57.51 

0.0177 

p-Cresol 

35.00 

62.32 

0.0163 

Pentadecane 

39.50 

70.33 

0 . 011 * 1 * 

Hexadecene 

1*6.00 

81.90 

0.012U 



RANGE OF THE DISPERSION PARAMETER Pe 
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0 

D 
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Column 

Carrier 

Peclet Number 

Dimensions 

Gas 

Range 

0.25mm dia. x 

helium. 

ll*6 - 161* 

0.025m long 

0.5ml /min 

, 

6.3mm dia. x 

helium. 

3560 - 1*910 

2.0m long 

60ml/min 

• 

0.76mm dia. x 

helium. 

16000 - 30000 

6.1m long 

2ral/min 


2.5mm dia. x 

helium, 

61*00 - 15000 

l*.9m long 

1* Oml /min 
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CARLE MINI -VALVE OPERATION 


BEFORE INJECTION 


IOOrI.LOOP 



CARRIER 6AS 


TO DETECTOR 


AFTER INJECTION 


IOO ill. LOOP 



CARRIER CAS 
TO DEI ECTOR 


Fig. 86. Concept of minivalve for injecting gas samples, 
















detailed in a recent report, Ref. 37, have been summarized below. 

Two versions of the equilibrium adsorption model were investigated. The 
first version, Eqn. 1 of Fig. 90 Incorporating a linear isotherm, has been 
solved analytically, Ref. 20, 21. The second version, Eqn. 2 of Fig. 90, 
incorporated the Langmuir isotherm into the model as a means for accounting 
for the observed composition-dependent effects. Analytical techniques were 
not available for solving the resulting nonlinear equation, so numerical 
methods were used. Because of instability problems with an orthogonal 
collocation procedure used previously, Ref. 34, the simple, finite differ- 
ence scheme of Crank-Nlcolsen was used. To obtain experience with the method, 
it was first applied to the linear version, Eqn. (1), since an analytical 
solution was available for comparison. The Crank-Nicolscn method replaced 
partial derivatives by finite differences, Fig. 91, which can be viewed as 
double averages of the compositions about the point of interest: position, 

(i); time, (N + 1/2). These finite difference expressions were substituted 
into the partial differential equation, giving a system of algebraic equations 
having a tridiagonal system matrix which was solved by the Thomas algorithm. 
Ref. 28. Considerable effort was made in the coding to eliminate calculations 
involving multiplication by ;:ero in order to reduce computation time. This 
additional coding reduced the computing effort by about 35%. Good comparison 
between the numerical and analytical results gave confidence in the method. 

To apply the method to the nonlinear version, Eqn. (2), a modification was 
required because of the nonlinear coefficient of the time derivative. A 
predictor-corrector technique was used. Initially the partial derivatives 
were replaced by finite differences as in the linear version. The composition- 
dependent coefficient of the time derivative was treated by a method developed 
by Douglas, Ref. 29. As shown in Fig. 92, yi >n +l/2 approximated by its 

value at the previous time step, yj n and a correction factor ( T/2) 

( y/ T) i>n where T is time. This derivative was approximated by its 
value from the partial differential equation in which the derivatives and 
compositions were evaluated at the previous time step. This approximation 
for the nonlinear term was substituted into the differential equation, 

Eqn. 2, and the resulting system of algebraic equations was solved by the 
Thomas algorithm. A new estimate for the nonlinear term was then made by 
averaging the compositions 


y l,n + 1/2 


(y l,n + y l ,ntl ) 12 


This result was substituted into the partial differential equations, and the 
procedure was continued until convergence was obtained. The technique pro- 
vided a workable means for numerically solving Eqn. 2 and evaluating it as a 
chromatograph model. 

The effect of changing adsorption parameters in this simulation model have 
been shown in Fig. 93 for three cases, each having the same peak time. The 
chromatogram for the simple, linear adsorption model (K ■ 0) was sharp and 
exhibited little spreading. As K Increased, the leading edge remained 
sharp while the trailing edge tailed. This type of spreading was typical 
of data not represented well by simple linear model. It was noted that the 
Langmuir constant K affected the peak time, so the thermodynamic parasKter 
oRq was reduced as K was increased to maintain a constant peak time. 
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The chromatograms of Fig. 93 were prepared to simulate data for n-heptane 
on Chrotnosorb C-102 obtained earlier, Ref. 34. To fit these data, a 
secant procedure adjusted the values of K to provide the observed peak 
time for an assumed ioRq. If predicted and observed gas compositions at 
the peak time differed, a new value of mR^ was selected, and the pro- 
cedure was repeated. Table 8 has summarized results for several simula- 
tions of the n-heptane data. Good representation of the data was obtained 
for valuesof 0.0426 and 8.6 for mRg and K respectively. Since the 
nonlinear model also reduced to the linear model for dilute solutions, it 
was concluded that it offered promise as a general simulation tool for the 
chemical systems encountered in this study. 

The Crank-Nicolsen procedure, while providing good numerical solutions 
to the nonlinear model, was computationally inefficient. Depending upon 
the chromatographic system studies, between 25 and 60 minutes of IBM 360/67 
computer time were required to fit a set of experimental data. Such com- 
putational efforts were excessive for practical application of the model, 
so future work will concern the development of a more efficient numerical 
procedure. 
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. \ TABLE 8 


COMPARISON OF THE SIMULATION PEAK HEIGHTS FOR 
VARIOUS MHO'S WITH THE ACTUAL PEAK HEIGHT 
(N-HEPTANE ON CHROMOSORB-102 AT 200°C) 


Thermodynamic 
Parameter, mRo 

Langmuir 
Constant, K 

Peak Height 
(mole fraction) 

0.0382 

12.2 

0.0308 

0.0426 

CO 

• 

cn 

0.0333 

0.0482 

*b. 

* 

oo 

0.0430 

ACTUAL DATA 


0.0353 
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THE EDUCATIONAL IMPACT OF NASA GRANT NGL 33-018-091 
AT RENSSELAER POLYTECHNIC INSTITUTE 


I. INTRODUCTION 

The project supported by NASA Grant NGL 33-018-091 was Initiated in 
the Fall of 1967 and has continued to the present time. During this period 
of time the project goals have been modified to reflect technical progress 
at Rensselaer and elsewhere and in reaction to long-range NASA planetary 
goals. During the past four years, the project goals have been consolidated 
and focussed on problems related to the "housing, protecting, servicing and 
transporting" of scientific instrumentation by an unmanned roving vehicle 
for the exploration of the planet Mars. 

As of this time, an imaginative vehicle concept has been conceived, 
designed, constructed and tested. In depth studies of terrain sensing con- 
cepts and devices, of obstacle (hazard) detection systems and of path 
selection systems simulation, development and appraisal, have been conducted 
and are continuing. All of these tasks focus on technical questions which 
NASA has judged to represent prime challenges to successful execution of aa 
unmanned exploration of the planet Mars. These tasks and the additional 
program in optimal design of a gas chromatograph separation system were 
judged to be germane not only to NASA but also to the team of faculty members 
involved. 

In appraising the value of this effort at Rensselaer over tie past 
eight years two of NASA's original goals in establishing the program should 
be considered. One of these goals was to obtain technical support in develop- 
ing a vehicle and related systems which would serve as a roving laboratory for 
unmanned exploration of Mars. The question of value with respect to this goal 
can be judged by individuals with technical competence and responsibilities. 

An indication of the productivity of the effort may be gained from the listing 
of reports and papers and presentations contained in Appendices A and B. 

This summary report addresses the second of NASA's goals, namely, to 
provide meaningful educational experiences in design for enginee* ing students 
at all educational levels. It is our conclusion that the program, which has 
Involved some 188 student participants at all academic levels, has been an 
outstanding success in this regard. The basis for reaching this conclusion 
and a description of how the objectives of engineering design education are 
achieved are described below. 


II. THE EDUCATIONAL PROCESS 


In establishing the specific technical tasks, three criteria were 
employed : 

1. The tusks must be relevant and critical to 
projected NASA missions, 

2. The tasks must offer significant opportunity 
for the professional growth of the faculty in- 
volved and must be in areas where the faculty 
have special expertise, 

3. The tasks must involve problems providing 
meaningful design experiences for participating 
students. 
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These criteria have guided the faculty team in proposing technical 
activities. During the past eight years some ten malor tasks, Table I, 
have been explored and a very much larger number of sub-tasks have been 
pursued. Some 188 different students have participated in these various 
activities for periods of time ranging from one semester to three or more 
years. In appraising the value of the educational experience, the organi- 
zation of the project, the student responsibilities, the reporting system, 
the role of the program monitors as well as the design process are worth 
reviewing. 


i 

II* A. Project: Crganizat ion 

The typical overall project organization is illustrated in Figure 1 
which applied to the 1973-1974 academic year* The five major tasks are 
coordinated by a project director who does not have direct responsibility 
for the work of individual students. However, the project director working 
with the faculty responsible for each major task area, and at appropriate 
times with student participants, defines the major scopes, goals, targets 
and inter-relationships. The faculty responsible for the major tasks work 
with the student to define more detailed scopes, goals and targets consistent 
with the overall plan. Figure 2, which is typical, illustrates the organiza- 
tion of the vehicle tasks. Sub- tasks are themselves sub-divided and organized 
according to the needs of the sub-task. 

One of the project directors major responsibilities is to identify 
’'interface" activities between major tasks and to develop those inter-relation- 
ships between sub-task groups required to insure compatability of the technical 
effort. As noted in a later section, this aspect contributes in a very real 
way to the student’s education in that the engineer is generally not free to 
optimize his solution to a particular problem without considering the impact 
of his proposal on other components of the system. 

This organization also emphasizes the need for the two-way flow of 
information to assure that tasks are proceeding on schedule towards a con- 
solidated goal. 

il.B, Student Responsibilities 

The student’s responsibilities include: definition of his problem 

within the framework of the major goals which have been specified, specifica- 
tion of constraints and performance requirements, development of a schedule 
for achieving the defined goals, performance of the indicated task within the 
specified time and preparation of appropriate reports. These responsibilities 
represent a radical change in the student’s academic experience. For the 
first time in his career, he is posed with defining his problem on the basis 
of reality and looking for imaginative and optimal solutions which are 
genet ally severely constrained. While it Is trie that he has the benefit of 
faculty review rmd guidance, he has to carry the main burden. Preparation of 
and adherence to a reasonable time schedule for a task which may last a full 
academic year or more is also a new and valuable experience as are the identi- 
fication of key barriers to an acceptable solution and the perception of problem- 
solving as a planned and logical sequence. Not the least of his new experience 
is the need to convey through oral and written reports the consequences and 
worthiness of his investigations. 
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The student finds that his academic background rarely provides all 
the knowledge and tools that he will need to solve his problem. He must 
gain additional knowledge from textbooks, handbooks, venoor's materials, and 
technical journals. He may have to conduct special research to gain the 
necessary information. He will have to build or find alternatives to 
demonstrate the validity of his recommendations. This is sharp contrast 
to his regular course work where his assignments are inevitably tailored 
to the subject matter. 

The student also finds that frequently the problem requirements ex- 
tend beyond his nominal disciplinary interest. For example, many of the 
mechanical engineering students working on the vehicle frequently find that 
electrical and electronic considerations are of equal importance to the 
mechanical considerations in resolving the problem at hand. 

II. C. The Feporting Responsibilities 

The reporting responsibilities of t lie students are illustrated in 
Figure 3 which is a typical program schedule. The project work refers to 
the schedule developed by the student, his faculty advisor and the project 
director as discussed under Section II. B. The student is responsible for 
preparing brief written summary reports for the technical monitor prior to 
each of five or six progress review meetings. The student must also make 
an oral presentation using such audio-visual materials and demonstrations 
as appropriate at these review meetings. He is required to prepare material 
for the mid-year and end-of-year progress reports summarizing project pro- 
gress for the interval. Finally, he is responsible for either as an 
individual or as a member of a team for writing detailed technical reports 
dcsciJ the investigations at key milestones. These technical reports 
receive wide distribution within NASA and other governmental agencies and 
interested industrial groups. 

The educational benefit derived by these reporting experiences also 
contrasts sharply *’ith normal academic programming. The student realizes 
for the first time that his documentation has the serious objective of con- 
veying meaningful information to readers who are not necessarily familiar 
with the specific problem and who can be expected to be professionally 
critical. The faculty and the project director make it abundantly clear that 
the reporting is more than a nominal chore, and that quality reporting is vital 
not only for this project's purposes but also as a routine requirement of 
quality professional practice as an engineer. 


I I . D. The Role of the Program Monitors and Formal Presentations 

During the past eight years, the students have benefited greatly from 
the involvement of program monitors representing the Jet Propulsion Laboratory. 
These monitors and other visitors selected on the basis of their particular 
expertise are present for the periodic progress review meetings at which the 
students make oral presentations of their work. The monitor has a multiple 
responsibility. He is a critic in the sense that he appraises the technical 
quality of the work. He also serves as a colleague in suggesting alternative 
approaches, methods, materials, etc. He provides valuable guidance In helping 
direct the work to useful purposes and is an important source of information. 


4 


From t.ho student's view, the monitor is a client. The student's 
goal is to transmit clearly and concisely the state of his work to the 
monitor and to be able to respond effectively to questions, criticisms 
and suggestions. 

The impact of the technical monitors, the selected guests, and the 
industritl ax.d governmental visitors at the annual review is profound. 

The student realizes that he is facing a "hostile" audience and is no 
longer in the sheltered academic womb. The audience is not personally 
"hostile" but it is an audience of practicing professionals whose personal 
careers and challenges demand a certain level and quality of performance. 

The student realizes at an early point that his work, is going to be judged 
on this more serious level. This realization has a substantial effect on 
how the student approaches his particular task and the intensity with which 
he applies himself as major milestones loom in the schedule. 

II. E. The Design Process 

Although there are many ways to describe the design process, the 
elements of the process inevitably include four major steps: 

1. Problem definition, constraint identification 
and information gathering. 

2. Development of alternative concepts. 

3. Analysis of alternatives. 

4. Decision making and implementation. 

The process is inherently iterative since new information and know- 
ledge may require a restatement of the problem, may offer new alternatives 
or may demand additional analysis. The objective is, of course, to try to 
minimize this activity but above all to move towards the strongest solution 
within the time constraints. 

The student's experience implants firmly this type of a design pro- 
cedure into his personality. He finds that these practical problems do not 
yield well to a haphazard attack jhe finds by experience and example that 
they can be approached systematically and creatively. His faculty advisor's 
questions, the project director's inquiries and the monitor's comments re- 
inforce dramatically the concepts of: Is this really the problem? Will 

your concept be compatible with the remaining system? Will It meet the 
established goal? Are there not better ways? etc. It is observed that the 
student's attitude and operational method have been strongly Influenced by 
participation in this project. 


III. STUDENT INVOLVEMENT AND PRODUCTIVITY 

It was noted earlier that the participating faculty have concluded that 
the NASA's educational objectives in supporting these activities have been 
achieved to an exceptional degree. It is fair to Inquire as to the basis of 
th:*s conclusion since it is in the self-serving Interest of the faculty in- 
volved to reach this favorable conclusion. 


One measure of the quality of the educational experience in design 
lies in the solutions to the problems as proposed in the technical reports 
and/or represented by the hardware and software produced. On an overall 
basis, the faculty appraisal is very favorable. Although a minority of 
the students do not achieve the desired educational goals, the majority do 
demonstrate a substantial growth and maturation and some develop to an ex- 
ceptional degree as evidenced by the technical reports listed in Appendix A. 
The faculty view appears to be supported by the favorable appraisals of the 
monitors and guests not only of the quality of the technical work but also 
of the presentations. 

An alternative measure of the quality of the technical work is pro- 
vided by the papers which have been published and the presentations which 
have been made at professional meetings. The diversity and extent of 
this activity is exhibited in Appendix B listing some AO papers and pre- 
sentations, the large majority of which are critically reviewed. This 
would suggest that the quality of the work is regarded as being professionally 
acceptable. It is worth noting that most of the papers are co-authored by 
faculty and students providing an additional enrichment of the students 
academic experience. 

The quality of the student's work is also reflected by the fact that 
Rensselaer students received Second Prizes in the Graduate Divisict of 
contests sponsored by the Lincoln Arc Welding Foundation in both 1970 and 
1971. 


Perhaps the most compelling argument of all is the enthusiasm of 
students participating in the project over the years a evidenced by the 
number involved. Table II lists the number of parcicipating students by 
their educational level for each of the eight years of the project. These 
data are also shown graphically in Figure 4. Some 188 students have had a 
significant involvement in the project including 8 doctoral degree, 110 
master degree and 70 baccalaureate recipients. The involvement ranged from 
a minimum of one semester (relatively rare) to three years for most doctoral 
candidates. The average involvement of master's and bachelor's candidates is 
about one full academic year. 

The strong student response, which has been growing as shown in Figure 
4, is due primarily to the student's perception of project participation as 
a valuable is well as interesting educational experience. The very large 
majority of the students participate without receiving any financial support 
from project funds since the project can only provide such assistance for 
5.5 to 6.0 full-time equivalent research assistants/yr or a maximum of 48 
such assistantships over the years. Since the doctoral candidates probably 
represent some 20 F.T.E. assistantships during their residence, it is i lear 
that the attractiveness of the project is due more to the atudent's perception 
of the worthiness of the experience and training and less to the prospects of 
financial support. 


IV. DISCUSSION 

In appraising the value of the project several questions deserve con- 
sideration: 


1. Has the attitude and ability of the participating 
students to address and resolve engineering design 
problems been developed to the desired level? 


2. Would these students have achieved these qualities 

in the absence of the project? \ 

3. What have been the professional careers of these 
students following graduation? 

A. What are the technological implications of this 
type of a NASA program? 

Faculty related to the project and other faculty aware of the activity 
are convinced that the attitudes and abilities of the students towards 
engineering design have been strengthened substantially and that the majority 
of students take with them on graduation a unique perspective towards pro- 
fessional practice. This conviction, which is based on our intimate inter- 
actions with the students, is supported by the contents of the technical 
reports, publications and presentations listed in Appendices A and B. It is 
also supported by the reactions of our program monitors, guests and visitors 
and by the awards our students have received. 

This experience in design provided to Rensselaer students is particularly 
germane since design education at engineering schools across the nation con- 
tinues to represent a low priority to the academic faculties on average. 
Despite the pressures of the accreditation societies to restore vigorous, 
Imaginative and in-depth education in design to its rightful place over the 
past fifteen years, strong programs in design remain relatively few in number 
and most design education is severely lacking in reality and relevance. This 
is not to infer that engineering students are not exposed to a reasonable 
number of courses that relate to design but rather that the emphasis is on 
standard approaches to well-defined areas within the disciplines involved. 

The students rarely have an opportunity to address open-ended problems with 
relevance and real need. It is our conviction that the large majority of 
the 188 students involved in this project would not have achieved these 
educational benefits in the absence of the NASA project. 

The dilemma of providing meaningful experiences in design has been pro- 
claimed in numerous publications over the past decade. The di lemma appears 
to be the result of two factors: (a) the real or imagined reward systems in 
engineering schools which emphasize scholarship in fundamental research areas 
and (b) the lack of real problems which are in fact relevant to some element 
of our society. These two factors arc not unrelated. The first factor can 
be made less significant or be eliminated depending on the extent that in- 
dustry and governmental agencies see fit to recognize and support relevant 
problem areas. Truly significant problems inevitably require faculty involve- 
ment of quality and this involvement will lead to new concepts and methods 
worthy of academic recognition. 

It is in this respect that NASA's generous support over the years has had 
a profound impact on education at Rensselaer. Faculty participants do not 
have to be coerced to become involved. Their labors are rewarded not only by 
the satisfaction of helping these young people mature but also more tangibly 
by the respect shown by colleagues and supervisors. The project has not only 
produced good problem solutions for NASA but has also led to papers and pre- 
sentations enhancing the faculty member's image at Rensselaer and elsewhere. 
The enthusiastic participation of faculty carries over to the students who 
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become attracted to the project. In Cum, this reaction has led to widespread 
recognition of the project by the student body on the campus. During the past 
four years, th. "Mars Project" has been one of the highlight presentations to 
prospective students, at Parent's Weekends and at alumni affairs. The publi- 
city the project received on campus has been instrumental ir encouraging 
improvements in a monber of design coutses. Many students other than those 
participating on the project have gained in a more direct way because the 
participating faculty carry over these design concepts and experiences in 
course;, which they are teaching. Thus there has been a direct impact on the 
teaching of design. 

While these peripheral benefits are welcome, the long-lived impact has 
been on the students and faculty who have participated in project activities. 
Their attitudes, procedures and abilities have been molded in an irreversible 
manner. The key to these educational benefits has been NASA's willingness 
to provide: (1) problems of substance whose solutions are truly needed, and 

(. > the financial support to make the problem solving process possible. NASA's 
commitment to the project has had an unquestioned effect on faculty and 
student interests. The Jet Propulsion Laboratory's commitment in providing 
program monitors has had a similar effect. Without these conmitmcnts it is 
not believed that the in-depth involvement of Rensselaer faculty and students 
could have been maintained. Indeed, with NASA's moral and financial support, 
the project activity has intensified in time with the result that this year 
37 students are participating in the program. 

One of NASA's and other governmental agencies' major concerns with 
supporting educational activities at universities is that students may be 
encouraged to follow programs of study leading to specializations incompatible 
with national needs or interests. In some cases, these agencies have found 
it p' -essary to maintain research programs of basic long-term value but 
relatively low priority to provide reasonable employment for the resulting 
graduates. This has certainly not been the case of this project at Rens- 
selaer. Sunsnarized in Table 111 are the employment data of degree recipients 
over the three-year period, July 1, 1971 through June 30, 1974. Of the 41 
graduates for whom we have data, all but three selected positions in the 
civilian economy. Only a small number of graduates were Involved with aero- 
space concerns reflecting the reduction of emphasis on the spec program in 
recent years. While we cannot document their activities with precision, 
the contacts we have had with a considerable number of these graduates indicate 
that the majority are using and benefiting from the skills they gained from 
the project. Samplings of the post graduation activities of graduates of 
earlier years follow similar patterns. Thus, this project aimed at producing 
technology for NASA har ’ ad the by-product effect on contributing to the 
civilian economy. 

Although the question of technologlcel implications of this type of 
NASA program is not directly germane to the educational question, it is 
worthy of consideration. The vast majority of research projects at univer- 
sities attract only those students for whom Che projects provide stipends 
and tuition. Aa noted earlier, two-thirdr and more of the effort has been 
provided by students who do not receive any direct financial support from 
the project. Their involvement is due to their Interest in the individual 
problems and the project as a whole and to what, we believe, is their per- 
ception of the educational value of participation. Thus a project of this 
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type has the effect of '^nultiplying" NASA funds. A relatively modest level 
of support produces student and faculty efforts far greater than what can 
be expected on a direct quid pro quo basis. The technical productivity in 
terms of hardware and software far exceeds that which could reasonably be 
expected from a normal "business" relationship. Thus, over the years NASA 
has gained a substantial "leverage" on its funding in terms of the producti- 
vity. 


The creativity inherent in young people and in faculty represents 
another aspect of the technological considerations. While it would be 
admitted that established aerospace vendors are without question far more 
qualified to design and produce the space-qualified hardware and software 
required for the missions, we would challenge with considerable vigor that 
such vendors can develop concepts more creatively. Indeed, we have some 
conviction that most industrial organizations cannot compete in this regard. 
As an example, the roving vehicle concept which has been developed over the 
past four years at Rensselaer is considered by many outsiders as a signifi- 
cant candidate for a planetary rover. It has some extraordinary features 
which are not possessed by any of the concepts proposed by industrial con- 
tributors. Neither the students nor the faculty are constrained by precon- 
ceptions as to what may be possible. The net effect is that the devices and 
systems which are proposed incorporate many novel features. 

In summary, the project has provided an exceptional experience and 
opportunity for students and faculty. Both groups have been enriched signi- 
ficantly by the activity. A significant number of unusually qualified 
graduates have entered professional practice in the civilian sector and 
substantial technical contributions have been made to NASA and the public. 


V. CONCLUSIONS 


1. The objective of developing the ability of engineering 
students to resolve design problems of substance has been 
achieved without qualification. The students who have 
participated on the project have gained an experience 
which is rare at universities. 

2. In the absence of the NASA support, these students 
would not have developed their ability in design to the 
levels which were achieved. 

3. The eight faculty members who have been involved directly 
have gained valuable educational experiences which have 
carried over in their teaching and have impacted on the 
education of many more students. 

4. The project activity attracted some 188 student parti- 
cipants including 8 doctoral, 110 masters and 70 bachelor 
degree recipients over the past 8 years. Less than one- 
third of the participating students received stipend 
and/or tuition support from the project. 

5. The project's success in attracting such large manber 
of students is due primarily to: 


(a) NASA’s commitment in providing ’’real 1 * problem 
opportunities and the support for seeking 

. further solutions, and 

(b) JPL’s commitment in providing interested and 
effective program monitors. 

The very large majority of graduates entered the civilian 
economy and are contributing to the national interest in 
widely diverse fields. 

The technological productivity of the project far ex- 
ceeded expectations for the funding levels involved. 

This is due to the large number of students participating 
without compensation. Thus, a very substantial ’’leverage* 1 
of NASA funds was obtained. 

The technical contributions to the NASA programs and to 
the nation at large were substantial and are reflected 
in hardware and roftware which was produced. These re- 
sults have been disseminated to the public through almost 
50 technical reports and some 40 papers and presentations 
at professional society meetings. 

The inherent creativity of young people is amply demon- 
strated by the hardware and software which was developed 
in response to the problems which were addressed. 
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Figure 3. Typical Program Schedule 











TABLE I 


SUMMARY OF MAJOR TASKS 


On-Line Updating of Martian Atmosphere Parameters during Entry 
Adaptive Trajectory Control 
Feasibility of Autogyro Landing System 
Global Navigation of a Mars Roving Vehicle 

Optimization of an Unmanned Mars Surface Exploration Mission 
Design, Construction and Evaluation of a Mars Roving Vehicle 
Obstacle Detection and Terrain Modeling for an Autonomous Rover 
Hazard Detection Systems 

Path Selection Systems Simulation, Development and Evaluation 
Chemical Analysis Systems Design 



is estimated. 


TABLE III m 

EMPLOYMENT OF PARTICIPANTS FOLLOWING GRADUATION, 1972- 1974 1 ' 


IBM . 5 

General Electric 5 

Raytheon 4 

Eastman Kodak 3 

Rohm and Haas 2 

Armed Forces 1 

U.S. Air Force Academy 1 

Construction Ind. 1 

Belfatcs & Pavarini, Arch. 1 

Turnkey Systems, Inc. 1 

Proctor and Gamble 1 

Burroughs 1 

Rochester Telephone 1 

Sanders Associates 1 

Union Carbide 1 

General Motors 1 

Pratt and Whitney 1 

Sikorsky 1 

Babcock & Wilcox 1 

Bell Laboratory 1 

Naval Architecture Lab 1 

Stone and Webster 1 

• Hewlitt Packard 1 

Corning Glass 1 

U.S. Government Lab 1 

Western Electric 1 

Returned to India 1 
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The total number of degrees earned by participants during the 
academic years, July 1, 1971 through June 30, 1974 la estimated 
at 57. Since alumni employment records are seldom complete, our 
knowledge of employment following graduation is limited to 41 of 
these students. 


APPENDIX A 


KPi TErKNTr*.^ REPORTS 


Sliva, Thomas F., "Chromatogiu^hir c *.-j fc em® s r s : First-Order Model Evaluation," 

RPI Technical Report MP-1, Rensselaer Poiytccnnic Institute, Troy, N.Y., 
September 1968. 

Jaoosko, R.E. , "On-Line Parameter Updating of the Martian Atmosphere with Minimum 
Storage," RPI Technical Report MP-2, Rensselaer Polytechnic Institute, Troy, 
N.Y. , January 1968. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and 
Surface Vehicle Control System for Mars Exploration," RPI Progress Report 
MP-3, Rensselaer Polytechnic Institute, Troy, N.Y., February 1969. 

Frederick, D.K. , et al. , "Analysis and Design of a Capsule Landing System and Sur- 
face Vehicle Control System for Mars Exploration," RPI Progress Report MP-4, 
Rensselaer Polytechnic Institute, Troy, N.Y., July 1969. 

Hedge, L.F., "Trajectory Control for Mars Entry by Discrete Changes of Drag Sur- 
face and Flight Path," RPI Technical Report MP-5, Rensselaer Polytechnic 
Institute, Troy, N.Y., August 1969. 

LaBarbera, G. , "Short Range Obstacle Detection Systems," RPI Progress Report MP-6, 
Rensselaer Polytechnic Institute, Troy, N.Y., September 1969. 

Rayfield, W.P. , "The Design Study of a Lightweight Rotor Hub for a Mai s- Landing 

Autogyro," RPI Technical Report MP-7, Rensselaer Polytechnic Institute, Troy, 
N.Y., September 1969. 

Voytus, W.A. , "Chromatographic Systems Analysis: Moment Analysis of the Equili- 

brium Adsorption Model," RPI Progress Report MP-9, Rensselaer Polytechnic 
Institute, Troy, N.Y. , August 1969. 

Krum, R.C., "Chromatographic Systems Analysis: Sample Injection Proolem," RPI 

Progress Report MP-10, Rensselaer Polytechnic Institute, Troy, N.Y., August 
1969. 

Mancini, R.J., "Terrain Modeling for Surface Vehicle Guidance," RPI Technical Re- 
port MP-11, Rensselaer Polytechnic Institute, Troy, N.Y., August 1969. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and 
Surface Vehicle Control System for Mars Exploration," RPI Technical Report 
MP-12, Rensselaer Polytechnic Institute, Troy, N.Y., March 1970. 

Wilson, J.V., "Instrumentation Study of Primary Navigation System for a Mars Roving 
Vehicle," RPI Technical Report MP-13, Rensselaer Polytechnic Institute, Troy, 
N.Y. , May 1970. 

Pavarini, C., and Chrysler, Jr., J.H., "Terrain Modeling and Path Selection by an 
Autonomous Martian Exploratory Vehicle," RPI Technical Report MP- 14, Rens- 
selaer Polytechnic Institute, Troy, N.Y., June 1970. 

Frederick, D.K. , et al., "Analysis and Design ot a Capsule Landing System and 
Surface Vehicle Control System for Mars Exploration," RPI Technical Report 
MP-15, Rensselaer polytechnic Institute, Tioy, N.Y., July 1970. 









i 1 - -J 



RPI TECHNICAT. REPORTS 'Continued) 


Rayfield, W.P., "Mars Roving Vehicle Conf iguralion, " RPI Technical Report MP-16, 
Rensselaer Polytechnic Institute, Troy, N.Y., August 1970. 

Goldberg, A.L. , and Frankel, P.J., "A Body-Bound Navigation System," RPI Technical 
Report MP-17, Rensselaer Polytechnic Institute, Troy, N.Y., June 1970. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and 
Surface Vehicle Control System for Mars Exploration," RPI Technical Report 
MP-18, Rensselaer Polytechnic Institute, Troy, N.Y., January 1971. 

Baer, S.R. , and Benoit, G.L. , "Chromatographic Test Facility," RPI Technical Report 
MP-19, Rensselaer Polytechnic Institute, Troy, N.Y., March 1971. 

Rieback, D.J., "Preliminary Design of an Automatic Device for the Location of the 
Pole Star and/or True Pole of Mars," RPI Technical Report MP-20, Rensselaer 
Polytechnic Institute, Troy, N.Y., May 1971. 

Simon, R.L. , "Design of a Toroidal Wheel for a Martian Roving Vehicle, RPI Techni- 
cal Report MP-21, Rensselaer Polytechnic Institute, Troy, N.Y., May 1971. 

Benoit, G.L. , "Reduction of Chromatographic Data and Evaluation of a GC Model," RPI 
Technical Report HP-22, Rensselaer Polytechnic Institute, Troy, N.Y., June 1971. 

Frederick, D.K. , ct al., "Analysis and Design of a Capsule Landing System and Sur- 
face Vehicle Control System for Mars Exploration," RPI Technical Report MP-23, 
Rensselaer Polytechnic Institute, Troy, N.Y., July 1971. 

Pavarini, C. , Baker, J. , and Goldberg, A., "An Optimal System Design Process," RPI 

Technical Report MP-24, Pensselaer Polytechnic Institute, Troy, N.Y, , Nov. 1971. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and Sur- 
face Vehicle Control System for Mars Exploration," RPI Technical Report MP-25, 
Rensselaer Polytechnic Institute, Troy, N.Y., July l to Dec. 31, 1971. 

Zuraski, C.D., "Laser Range Measurement for a Satellite Navigation S:heme and Mid- 
Range Path Selection and Obstacle Avoidance," RPI Technical Report MP-26, 
Rensselaer Polytechnic Institute, Troy, N.Y., June 1972. 

Keba, P.S., and Woodrow, P.T. , "Comparison of Two Gas Chromatograph Models and 

Analysis of Binary Data," RPI Technical Report MP-27, Rensselaer Polytechnic 
Institute, Troy, N.Y., July 1972. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and Sur- 
face Vehicle Control System for Mars Exploration," RPI Technical Report MP-28, 
Rensselaer Polytechnic Institute, Troy, N.Y., July 1, 1971 to June 30, 1972. 

Bohelm, S.L. , and Pur don, W.C., "Path Selection System Simulation and Evaluation for 
a Martian Roving Vehicle," RPI Technical Report MP-29, Rensselaer Polytechnic 
Institute, Troy, N.Y., December 1972. 

Frederick, et al., "Analysis and Design of a Capsule Landing System and Surface 
Vehicle Control System for Mars Exploration," RPI Technical Report MP-30, 
Rensselaer Polytechnic Institute, Troy, N.Y. , 3uly 1, 1972 to Dec. 31, 1972. 


\ 

i 




I 


I 

, i 



I 

i 


3 , 


RPi TECHNICAL REPORTS (Contittud) 


Heisch, A,J., "Binary Chromatographic Data and d^tissat *,r». of i'Jso.-’-.-nt Porosities," 
RPI Technical Report MP-31, Rensselaer Polytechnic Institute, Troy, N.Y., 

July 1972. 

Pfeifer, W.J., "Gradient Estimates from Stereo Measurements for a Martian Vehicle,” 
RPI Technical Report MP-32, Rensselaer Polytechnic Institute, Troy, N.Y., 

May 1973. 

Burger, Paul A., "Stochastic Estimates of Gradient from Laser Measurements for an 
Autonomous Martian Roving Vehicle," RPI Technical Report MP-33, Rensselaer 
Polytechnic Institute, Troy, N.Y., May 1973. 

Herb, G.T. , "Laser Scanning Methods and a Phase Comoarison, Modulated Laser Range 
Finder for Terrain Sensing on a Mars Roving Vehicle," RPI Technical Report 
MP-34, Rensselaer Polytechnic Institute, Troy, N.Y., May 1973. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and 
Surface Vehicle Control System for Mars Exploration," RPI Technical Report 
MP-35, Rensselaer Polytechnic Institute, July 1, 1972 to June 30, 1973. 

Woodrow, P.T., "Preliminary Numerical Analysis of Improved Gas Chromatograph Model," 
RPI Technical Report MP-36, Rensselaer Polytechnic Institute, Troy, N.Y., 

Sept. 1973. 

Palunbo, D.L. , "Design of a Laser Rangefinder for Martian Terrain Measurements," 

RPI Technical Report MP-37, Rensselaer Polytechnic Institute, Troy, N.7., 

May 1973. 

Ryder, Alan G. , "Dynamic Evaluation of RPI's 0.4 Scale Unmanned Martian Roving 

Vehicle Model, RP* Technical Report MP-38, Rf.ssclacr Polytechnic Institute, 
December 1973. 

Frederick, D.K. , et al., "Analysis and Design of a Capsule Landing System and 

Surface Vehicle Control System for Mars Exploration," Rensselaer Polytechnic 
Institute, RPI Technical Report MP-39, February 1974. 

Woodrow, P.T. , "Analysis of Chromatograph Systems Using Orthogonal Collocation," 

RPI Technical Report HP- 40, Rensselaer Polytechnic Institute, Troy, N.Y., 

May 1974. 

Lavoie, R.C. , "Composition Dependent Effects in Gas Chromatography," RPI Technical 
Report MP-41, Rensselaer Polytechnic Institute, Troy, N.Y., May 1974. 

Campbell, R.S. and Slmouds, R.R. , "Path Selection System Development anP Evaluation 
for a Martian Roving Vehicle," RPI Technical Report MP-42, Rensselaer Poly* 
technic Institute, Troy, N.Y., May 1974. 

Pavarinl, Carl, "System Design Optimisation for a Mars-Roving Vehicle end Perturbed- 
Optimal Solutions in Nonlinear Programming," RPI Technical Report MP-43, 
Rensselaer Polytechnic Institute, Troy, N.Y., June U74. 

Frederick, D.K. , et si., "Anslysis and Desig.'. of a Capsule Landing Sybtcm and 

Surface Vehicle Control System for Mars Exploration," Rensselaer Polytechnic 
Institute, RPI Technical Report MP-44, July 1, 1973 to June 30, 1974. 






0 

0 

u 


Reed, Katei:* A., "Kecier ition of Three Dimensional Obstacles by an Edge Detection 
Scheme, RP1 Technical Report MF-45, Rensselaer Polytechnic Institute, Troy, 
N.Y., May 1974. 

D'Angelo, Kenneth R. , "Parameter Estimation for Terrain Modeling from Gradient 
Data," RPI Technical Report MP-46, Rensselaer Polytechnic Institute, Troy, 
N.Y. , May 1974. 



L 

l 


L 

L; 

y 

D 

D 

8 

I 



4 


APPENDIX B 


RPI Pap;.* and Presentations 

f l Shen, C.N., and P, J. Cefola, "Adaptive Trajectory Control for Kars Entry B'sed on 
4 * Sensitivity Analysis," AIM 68-835, AIM Cuidance, Control and Flight Dy stales 
Conference, Pasadena, Calif,, August 1968. (Accepted for publication *n Che 
JV AIM Journal), 

Cefola, P.J. and Shen, C.N. , "Vector-Matrix Second Order Sensitivity Equation with 
jj Application to Mars Entry," AIM Journal, 7., No. 8, 1633 (1969). 

Shea, C.N. , and Cefola, P.J., "Adaptive Control for M'rs Entry* Based on Sensitivity 
, Analysis," AIM Journal, No. 6, 1145 (1969). 

‘ Shea, C.N. , and Janosko, R. , "On-Line Parameter Upo**-<ng of the Martian Atmosphere 

with Minimal Storage," Proceedings of the Second ilawaii International Conference 
on System Sciences at the University of Hawaii, January 1969. 

I i 

Broniarek, C.A. and Sandor, G. N., "Nonlinear Vibration Analysis of Four-Wheel 
j Surface Vehicle with Hydraulic Driving System," presented at and published in 

' the P roceedings of the Second International Congress on the Theory of Machines 

and Mechanisms, Zakopane, Poland, Sept. 24-27, 19.9, Vol. 2, pp. 45-55. 

. , Shen, A.N., and Yong, K. , "On-Line Tracking of Space Vehicles Usinj Modified Sequ:- 
* tial Estimation," Proceedings 3rd Hawaii International Conference on System 
Science, 1970. 

Sandor, G.N., and Rayfleld, W.P., "Dragster on Mars," ASME News, Hudson Mohawk 
Section, Vol. 7, No. 3, January- February, 1970, pp. 3-7. 

* Janosko, R.E. , and Shen, C.N. , "Consecutive Tracking of Landmarks: by an Active 

Satellite and then by a Martian P.oving Vehicle," Proceedings 3rd Hawaii Inter- 
| national Conference on System Science, 1970. 

Duke, J.C., and Shen, C.N. , "Optimal Filtering for Konlinca Parameter Updating of 
j the Martian A'mocpherc," Proceedings of the Symposium tor International Federation 

L Automatic Control, Kyoto, Japan, August 1970. 

I * Rayfleld, W.P. , and Sandor, G.N. , "Rcnssc lae-'s Roving Vehicle for Furs," presentee, 
at the First Western Space Congress, Santa Maria, Calif., October 27-29, 1970, 
published i* Proceeding, pp. 838-855. 

| Chen, H.M. , end Shen, C.K., "Surface Navigation System and Error Analysis for Martian 
** Rover," Proceedings of the Fourth Hawaii International Conference on Systra 
Sciences, January 1971. . 

I 

• Rayfleld, W.P., and Sandor, G.N., '‘Design of e Roving Vehicle for Mars," Invited 
Paper, presented at the ASME Design Engineering Conf., New York, N.Y. , 

| April 19-22, 1971. 

Cuke, J.S., and Shen, C.N., "On-Line Parameter Updating for Mars St.&to-.'phere and 

I Troposphere," Proceedings of the AIAA- Joint Automatic Control Conference at 
St. Louis, Mo., August 1971. 


I Cefola, P.J., and Shen, C.N,, "A Parameter Dependent Spacecraft Cuidance Boundary 
Value Problem," AIM Journal, Vol, 9, if 10, pp 1975-.1979, October 197 a. 

I Sandor, C.N., "Rensselaer's Roving Vtn^cle for Mars: lessons for Possible Cl*nlm*). 

Applications," Medical Society cC N.Y., 16th Annual Convention, .eb. 19/2, 




I 


l 




Jf . RPI Papers e-i ’’vesentatlons (continued) 

** Sandor, G.N,, "Hotorizeu Wheel Chair Can Even Cl < - K Sr-i.s," (German), Medical 
T ribune, an international weekly, German laneuaee issue roi West Germany, 
n Voi. 7, No, 11, Fri., Mar. 17, 1972, puhl<.,hcd in Wiesbaden, pp. 1 & 48; and 

j; also in the German language issue of th'» »eekl> p^oiisned in Baden for 

Switzerland, Vol. 5, No. 11, Fri., Mar. 17, 1972, pp. 1 '< 20. 

I ’ 

j; Sandor, G.N., "Wheel Chair Can Elevate to Eye Level," Chronic Disease Management 
Vol. 6, No. 4, April 1972, p. 18. 

J 1 Yerazunls, S., Moore, J.W., and Wehe, R. L. , "Unmanned Exploration of Mars: A 

• ■ NASA-Univcrsity Experiment in Engineering Education," American Chemical Society 

Industrial-Academic Cooperation Symposium, Boston, April 1972 (with Cornell 
and JPL). 

' ■ ' | 

Janosko, R.E., and Shen, C.N., "A Simplified Satellite Navigation System for an 
Autonomous Mars Roving Vehicle," Proceedings of the Fifth Congress of Inter- 
national Federation of Automatic Control at Paris, France, June 1972. 


t 


1 

L 

Li 

0 

11 

L 

I 


Broniarek, C.A., and Sandor, G.N. , "On Surface Vehicle Body Motion with Gyroscopic 
Storage of Kinetic Energy," Mechanism and * , achine Theory . Vol. 7, No, 2, 

Summer 1972, pp. 141-154. 

Pavarini, C., Smith, E.J., Van Donburg, N. , -’System Modeling and Optimal Design of 
, a Mars Roving Vehicle," Proceedings of IEEE Conference on Decision and Control, 

New Orleans, La., December 1972, 

Shen, C.N., and Burger, P. , "Stochastic Estimates of Gradients from Laser Measure- 
ments for an Autonomous Martian Roving Vehicle," Proceedings of 3rd IFAC 
Symposium, The Hague, June 1973. 

Sandor, G.N. , "Unmanned Martian Roving Vehicle Inspires Stairclimbing Wheel Chair," 

New Yo~k State Journal of Medicine, Dec. 15, 1973, Vol. 73, No. 24, pp. 2880-2883. 

Shen, C.N. , and D'Angelo, R.P., "Parameter Estimation for Terrain Modeling froa 

Gradient Data," 7th Hawaii International Conference on System Sciences, Hawaii, 
January 1974. 

Burr, A. A. and Sandor, G.N., "Teaching Problem-Solving at Rensselaer," invited paper, 
Journa l of Engineeri n g Education , February 1974, Vol. 64, No. 5, pp. 357-359. 
(Describes RPI's NASA Mars Rover Design Project). 

Reed, M, , Sanyal P, , and Shen, C.N., "A Practical Obstacle Detection System for the 
Mars Rover," Proceedings of the Milwaukee Symposium on Automatic Controls, 

March 1974, Milwaukee, Wisconsin, 

Sandor, G.N,, "Seven Dangers of Designer Overspecializaticn and How to Avoid Them 

by Designer Education," Invited ASME Paper No. 74-DE-35, 1974 Design Engineering 
Conference, Chicago, April 1-4, 1974, to be preprinted by the ASME in full and 
to be published in abstract in M echanical Engineering . Journal of the ASME. 
(Educational Atp«cts of RPI's NASA Mars hover Design Piojcct as Conducted in EDL - 
(RPI's Engineering Design Laboratory, Geof^e N, Sandov, Director. )) 

Sandor, G. N., "Seven Dangers of Designer Overspecialization and How to Avoid Them by 
Designer Education," Electromechanical Design, March 1974. 

Lashmot, P.K. and Woodrow, P.T., "Simulation of Pulsed, Distributed Systems Using 

\ Orthogonal Collocation," Salt Lake City Meeting, American Institute of Chemical 
Engineers, August 19 1 4. 

Lashmet, P.K., Benoit, G.L. and Woodrow, P.T. , "Evaluation of a Gas Chromatograph Model," 
Atlantic City Meeting, American Chemical Society, September 197a. 


3 . 


RPI Papers and Presentations (Continued) 


Sandor, G.N., Hassel, David R. and Marino, Philip F, , "Modern Mechanisms Make 
Manless Martian Mission Mobile — "Spiu-Off" Spells Stairclimbing Self- 
Sufficiency for Larthbound Handicapped", Proceedings of 9th Aerospace 
Mechanisms Symposium, 18-1 to 18-17, John F.* Kennedy Space Center, October 
1974. 

Shen, C. N. and J. S. Sher, "Simulation of the Range Data and The Picture Enhancing 
Scheme for the Mars Rover", ACM 1975 Computer Science Conference, Washington, 

D.C., February 1975. 

Shen, C.N. and Audbur Thompson, "Computer Results of Two Dimensional Spline Function 
in Terrain Modeling Optimization", ACM 1975 Computer Science Conference, 

Washington, D. C. , February. 1975. 

Shen, C.N. and Mark Friedman, "The Effect of Laser Measurement Reflection Error on 

Range Data Processing", Milwaukee Symposium on Automatic Computation and Control, 
Milwaukee, Wisconsin, April 1975. 

Shen, C.N. and R. V. Sonalkar, "Mars Obstacle Detection by Rapid Estimation Scheme 

From Noisy Laser Rangefinder Readings", Milwaukee Symposium or, Automatic Computa- 
tion and Control, Milwaukee, Wisconsin, April 1975. 

Shea, C.N. and M. Friedman, "Measurement Scanning Schemes for Terrain Modeling", 

6th Annual Pittsburgh Conference on Modeling and Simulation, Pittsburgh, Pa., 

April 1975. 

Shen, C. N. , Leung, K.L, Yerazunis, S. , "Classification of Terrain Models for Martian 
Vehicle," 6th Annual Pittsburgh Conference on Modeling and Simulation, Pittsburgh, 
Pa., April 1975. 

Frederick, D.K. , "Path Selection System Simulation and Evaluation for a Martian Roving 
Vehicle", 6th Annual Pittsburgh Conference on Modeling and Simulation, Pittsburgh, 
Pa., April 1975. 

Woodrow, P.T. and Lashmet, P.K, , "A Formulation Consideration for Orthogonal Collocation 
Procedures", AICA International Symposium on Computer Methods for Partial Differ- 
ential Equations, Bethelehem, Pa., June 1975. 

Woodrow, P.T. and Lashmet, P.K., "Zero-Derivative Boundary Conditions for Pulsed 
Distributed Systems", AICA International Symposium on Computer Methods for 
Parial Differential Equations, Bethlehem, Pa., June 1975. 


